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1 Introduction. 



In this paper we study triangular arrays of Markov chains X^h (0 < k < n) that 
converge weakly to a diffusion process Y s (0 < s < T) for n — > oo. Here h = T/n 
denotes the discretization step. We allow that T depends on n. In particular, we 
consider the case that T — > for n — > oo. Our main result will give Edgeworth type 
expansions for the transition densities. The order of the expansions is o(/i _1_l5 ), 5 > 0. 
This is done for time horizons that may converge to 0. The research is motivated 
by applications to high frequency data that are available on a very fine grid but are 
approximated by a diffusion model on a more rough grid. The work of this paper 
generalizes the results in Konakov and Mammen (2005) in two directions. The time 
horizon is allowed to converge to and also cases are treated with nonhomogenous 
diffusion limits. 

The theory of Edgeworth expansions is well developed for sums of independent 
random variables. For more general models approaches have been used where the 
expansion is reduced to models with sums of independent random variables. This is 
also the basic idea of our approach. We will make use of the parametrix method. In 
this approach the transition density is represented as a nested sum of functionals of 
densities of sums of independent variables. Plugging Edgeworth expansions into this 
representation will result in an expansion for the transition density. 

Weak convergence of the distribution of scaled discrete time Markov processes to dif- 
fusions has been extensively studied in the literature [see Skorohod (1965) and Stroock 
and Varadhan (1979)]. Local limit theorems for Markov chains were given in Konakov 
and Molchanov (1984) and Konakov and Mammen (2000,2001,2002). In Konakov and 
Mammen (2000) it was shown that the transition density of a Markov chain converges 
with rate 0(n -1 / 2 ) to the transition density in the diffusion model. For the proof 
there an analytical approach was chosen that made essential use of the parametrix 
method. This method permits to obtain tractable representations of transition den- 
sities of diffusions that are based on Gaussian densities, see Lemma 1 below. Similar 
representations hold for discrete time Markov chains X kth , see Lemma 3 below. For a 
short exposition of the parametrix method, see Section 3 and Konakov and Mammen 
(2000). The parametrix method for Markov chains developped in Konakov and Mam- 
men (2000) is exposed in Section 4. Applications to Markov random walks are given in 
Konakov and Mammen (2001). In Konakov and Mammen (2002) the approach is used 
to give Edgeworth-type expansions for Euler schemes for differential equations. Re- 
lated treatments of Euler schemes can be found in [Bally and Talay (1996a, b), Protter 
and Talay (1997), Jacod and Protter (1998), Jacod (2004), Jacod, Kurtz, Meleard and 
Protter (2005) and Guyon (2006)]. Standard references for the parametrix method are 
the books by Friedman (1964) and Ladyzenskaja, Solonnikov and Ural'ceva (1968) on 
parabolic PDEs [see also McKean and Singer (1967)]. 
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The paper is organized as follows. In the next section we will present our model for 
the Markov chain and state our main result that gives an Edgeworth-type expansion 
for Markov chains. In Section 3 we will give a short introduction into the parametrix 
method for diffusions. In Section 4 we will recall the parametrix approach developed 
in Konakov and Mammen (2000) for Markov chains. Technical discussions, auxiliary 
results and proofs are given in Sections 5-7. 

2 Results. 

Let n > 2 , T = T{n) < 1 and h = T/n. Suppose that q(t,x, •) , {t,x) E [0, 1] x R d 
is a given family of densities on R d , Xu(t,x) is the v — th cumulant corresponding to 
the density q(t, x, ■) and m is a function from [0, 1] x R d into R d . We shall impose the 
following conditions 

(Al) j Rd yq{t,x,y)dy = 0, < t < 1, x E R d . 

(A2) There exists positive constants cr* and a* such that the covariance matrix a (t, x) = 
fa* yy T( l (*' x iV) d V satisfies 

a* < 6 T a (t,x)9< a*, 

for all ||0|| = 1 and t E [0, 1], x E R d . 

(A3) There exists a positive integer S' and a real nonnegative function ip (y) , y E M. d 
satisfying sup^ggd ip (y) < oo and J Rd \\y\\ if) (y) dy < oo with S = (S' + 2)d + 4 
such that 

\r%q(t,x,y)\ <4>(y), te [0,1], x,y ER d \u\ = 0, 1, 2, 3, 4 
\D»q{t,x,y)\ < if){y), te [0,1], x, y E R d \v\ =0,1,2. 

It follows from (A2), Lemma 5 and (fTH| . (|T7]l below (applied with h = 1 ) that the 
following condition holds 

(A3') For all x, y E R, h > 0, < t, t + jh < 1, j > jo, with j does not depending 
on x, t 

\D v x q^ (t,x,y)\ < Cj- d/2 ^ {r 1/2 y) , W\ = 0, 1,2,3 

for a constant C < oo. Here q^'(t,x,y) denotes the j— fold convolution of q 
for fixed x as a function of y. 

q b) (t,x,y) = J q (j - 1 \t,x,u)q(t+ (j - l)h, X, y - u)du, 
q (1) (t,x,y) = q(t,x,y). 
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It follows also from (A3) that for 1 < j < jo 

\y\\ S q U \t 1 x,y)dy < C(j ). 



(Bl) The functions m(t,x) and a(t,x) and their first and second derivatives w.r.t. 
t and their derivatives up to the order six w.r.t. x are continuous and bounded 
uniformly in t and x. All these functions are Lipschitz continuous with respect 
to x with a Lipschitz constant that does not depend on t. The function Xv{t,%), 
\u\ = 3, 4, is Lipschitz continuous with respect to t with a Lipschitz constant that 
does not depend on x. A sufficient condition for this is the following inequality 

(1 + \\z\\ A ) \q{t,x,z) -q{t',x } z)\dz < C \t - t'\ , < t, t' < 1, 

with a constant that does not depend on x G M. d . Furthemore, D u £ a (t, x) exists 
for \u\ = 6 and is Holder continuous w.r.t. x with positive exponent and a 
constant that does not depend on t. 

(B2) There exist x < | such that liminf n ^ 0C T(n)n ,,r > (remind that we consider 
T(n) < 1 ). 

Consider a family of Markov processes in M. d of the following form 

(1) X k +i,h = X kjh + m (kh, X k ,h) h + Vh£, k+hh , X ,h = x G R d , k = 0,...,n-l, 

where (^i,h)j =1 n is an innovation sequence satisfying the Markov assumption: the con- 
ditional distribution of £fc+i,/i given X^h — x k, ■■■,X j l = x depends only on X k ^ = x k 
and has conditional density q (kh, x k , •) • The conditional covariance matrix correspond- 
ing to this density is a(kh,x k ). The transition densities of (X it h) i=1 n are denoted by 
p h (0,kh,x, ■) . 

We will consider the process (JIJ) as an approximation to the following stochastic 
differential equation in M. d : 

dY s = m (s, Y s ) ds + A (s, Y s ) dW s , Y = x(ER d , s G [0, T], 

where (W s ) s>0 is the standard Wiener process and A is a symmetic positive definite 
dx d matrix such that A (s, y) A (s, y) T = a (s, y) . The conditional density of Y t , given 
Y = x is denoted by p (0, t, x, •) . We shall consider the following differential operators 
L and L : 

Lf(st r v)--^a(s x) 9 ^^ *' X ' v) I Vm (o x) d f( s ^ x ^) 
Lj{s, t,x,y) - - o%j\.s,x) q x ,q x , + 2_^ m ^ s i x > dx . 

i,j=l 1 ^ i=l * 

(2) Lf(s,t,x,y) = l - <*ij(s,y) d f g x ,t d ^ y ^ +^ m i(^y) df ^ dx ,X,y " > - 

i,j=l 1 i i=l * 



4 



To formulate our main result we need also the following operators 

L'jis t x y) = 1 V s dai i( s > x ) d2 f( s ^^ x ^) + dmj(s,x) df (s,t,x,y) 
2 ^ ds dxidxj ds dxi 

i,3=l J 1=1 

/„> j, f ( . ,±^da ij {s,y)d 2 f{s,t,x,y) J^dm l (s,y)df{s,t,x,y) 
(3) L f [s, t,v,z)--^ — s + ^ — ^ ^ • 

and the convolution type binary operation (8> : 



/ ® g(s,t,x,y) = du f (s,u,x, z) g (u,t, z,y) dz. 

Js JR d 

Konakov and Mammen (2000) obtained a nonuniform rate of convergence for the 
difference ph (0, T, x, •) — p(0, T, x, •) as n — >• oo in the case T x 1. Konakov (2006) 
proved an analogous result for the case T = o (1) . Edgeworth type expansions for 
the case T x 1 and homogenous diffusions were obtained in Konakov and Mammen 
(2005). The goal of the present paper is to obtain the Edgeworth type expansions for 
nonhomegenious case and for both cases T x 1 or T = o(l) . The following theorem 
contains our main result. It gives Edgeworth type expansions for ph- For the statement 
of the theorem we introduce the following differential operators 



F x [f]{s,t,x,y) = ^*!^I%f(8,t,x,y), 



F 2 [f)(s,t,x,y) = £*^££/( a ,f,x, y ). 

[i/|=4 



Furthermore, we introduce two terms corresponding to the classical Edgeworth expan- 
sion 

(4) H^t,x,y) = (t-s)J2 Xu{S ^ y) D u x p(s,t,x,y), 

["1=3 



(5) Hs,t,x,y) = (t-s)J2 ^ \ ' U, Bffi8,t,x,v) 

["1=4 

+l(t-sf{j2^^D») p(s,t,x,y), 

where 



M=3 



l r l 

Xv(s,t,y) = - / Xu(u,y)du, 

and where Xv(t,x) is the v — th cumulant of the density of the innovations q(t,x, •). 
The gaussian transition densities p(s,t,x,y) are defined in (0). Note, that in the 
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homogenous case Xu(u,y) = Xu(y) and Xu( s ,t,y) = xAu), where xAu) is the u - th 
cumulant of the density q(y , ■). 

Theorem 1. Assume (Al)-(A3), (B1),(B2). Then there exists a constants > 
such that the following expansion holds: 



sup 

x,y£R d 



T d/2 1 + 



y-x 



T 



x |Pft(0; T, x, y) — p(0, T, x, y) 



-/^(O, T, x, y) - hn 2 (0, T, x, y) |] = 0(/i 1+5 ), 
where S' is defined in Assumption (A3) and where 

ni(0,T,x,y) = p®J 7 i\p])(0,T,x,y), 

7r 2 (0,T,x,y) = (p® F 2 \p})(Q,T,x,y)+p®F 1 \p® ^i[p]](0,T,a;,y) 

+ 7^ ® - £ 2 M0, T, x, y) - ~p ® (L' - L')p(0, T, x, y). 

Here p(s,t,x,y) is a transition density of the diffusion Y(t) and the operator is 
defined analogously to L but with the coefficients 'frozen" at the point x. The norm ||«|| 
is the usual Euclidean norm. 
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Discussion and remarks 

1. It can be easily shown that 

\h l/2 7n(0,T,x,y)\ < Cxn-V 2 T- d/2 exp 

\hir 2 {0,T,x,y)\ < C x rr l T- d ' 2 exp 
and that by definition of h 

0(h 1+s )T- 

with some positive constants C\ and C 2 . 

2. If the innovation density q(t, x, •) does not depend on x then — L,L' — V and 
p(s, t, x, y) = p(s, t, x, y) where p is defined in (JHJ) with a(s, t, y) = a(s, t) = J s a{u)du 
and m(s,t,y) = m(s,t) = J m(u)du. This gives 

7Ti(0,T,z,y) = I ds ! p{0,s,x,v) ^^-D"p(s,T,v,y)dv 

J J I I Q 









s'- 


-i 


d/2 


1 + 


y-x 










VT 













S'~ 


l-Srpl+S-d/2 


1 + 


y-x 






VT 





M=3 



E 

["1=3 



dsDy / p(0, s, x, v)p(s, T, v, y)dv 



T 



^X„(0,T)D»p(0,T,x,y) = v 1 (0,T,x,y), 



["1=3 



P<8> Fi\p\(s,T, z,y) = / du p(s,u,z,w) 

Js J \u\=3 V - 

= -J2 D v f T ' *> y) = ( T - s )J2 ^p-ms, T, z, y), 

I I n J S I I o 



|H=3 " a M=3 



^[p^^T^y) = (T- S ) ^ *j^D 



|i/|=3 



M=3 



= (t- S ) e ^^ rp(^,^), 

|H=3>'|=3 

p <g> ^ 2 [p] (0, T, x, y) + p ® ^ [p ® T x [p\ (0, T, x, y ) = T E ^pW (0, T, x, y) 
+ /^ds j p(0, W )(T- S ) E ^p^f^ fp^T,^) 

" / ° J M=3,M=3 



1 1,|=4 



Z/! 



M=3,M=3 

For v = v 1 we have 



So Xv<yS ^Ss Xu '^ du ) ds = \j Q J Q Xu(s)Xu(u)dsdu = yXv(Q,T)Xv(0,T). 
For v 7^ v 1 we consider 



n[Xi/(s)x„'(u) + X^'( S )X^(^)] rfs^M 

= \ J J ^ Xu ^ Xu '^ + xA s )Xv(u)] dsdu 
= yx„(0, ^(0, T) + yx„(0, T)x„(0, T). 
From the last equations we obtain 

P ® ^2 [p] (0, T, x, y) + p ® ^ [p ® ^ [p] (0, T, x, y) 



T 2 



X,(0,T) „ 



|H=4 

7r 2 (0,T,x,y). 



| „|=3 
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Thus for this case we get the first two terms of the classical Edgeworth expansion 
h^n^T^y) + hn 2 (0,T,x,y) for the sums of independent non identically dis- 
tributed random vectors. 

3. If Xu(t,x) = for H = 3 and for t G [0,T] x R d then it holds that J= x = 0. 
The Theorem 1 holds with 

7ri(0,T,x,y) = 0, 

n 2 (0,T,x,y) = (p <g) J r 2 [p])(0, T, x, y) 

+ l -p <g> {Ll - L 2 )p(0, T, x, y)- l -p® (L' - Z>(0, T, x, y). 

If in addition Xu(t,x) = for \v\ = 4 then the first four moments of the innovations 
coincide with the first four moments of a normal distribution with zero mean and 
covariance matrix a(t,x). In this case we have Ti = and 

7Ti(0,T,x,y) = 0, 

ir 2 (0,T,x,y) = X -p ® (L\ - L 2 )p(0, T, x, y) - l -p ® {U - L')p(0, T, x, y) 

and the first two terms of the Edgeworth expansion do not depend on the innovation 
density. In particular the Edgeworth expansion for the homogeneous Euler scheme 
holds with the same 7Ti and 7r 2 as in the last two equations. For the homogenous case 

ni(0,T,x,y) = 0, 

7r 2 (0,T,x,y) = ^p®(Ll-L 2 )p(0,T,x,y). 

This result for T = [0,1] under a weaker condition on the diffusion matrix was obtained 
by Bally and Talay (1996). 



3 Parametrix method for diffusions. 



For any s e [0, T], x, y e M, d we consider an additional family of "frozen" diffusion 
processes 

dY t = m (t, y) dt + A (t, y) dW t , Y s = x, s<t<T. 

Let p y (s,t,x, •) be the conditional density of Y t , given Y s = x. In the sequel for any z 
we shall denote p (s, t, x, z) = p z (s, t, x, z) , where the variable z acts here twise: as the 
arument of the density and as defining quantity of the process Y t . 

The transition densities p can be computed explicitly 



(6) 



p(s,t,x,y) = (2tt) d/2 (det a (s,t,y)) 1/2 

x exp ^-i (y - x - m (s, t, y)) T a' 1 (s, t,y)(y - x - m (s, t, y))^j , 



where 



a(s,t,y)= / a(u,y)du, m(s,t,y)= / m(u,y)du. 

J s J s 



Note that the differential operators L and L corresponds to the infinitesimal operators 
of Y or of the frozen process Y s ^ y , respectively, i.e. 

Lf(s,t,x,y) = \imh- 1 {E[f(s,t,Y(s + h),y)\Y(s)=x]-f(s,t,x,y)}, 



Lf(s,t,x,y) = \imh l {E[f(s,t,Y s , x , y (s + h),y)} - f(s,t,x,y)}. 

h— »0 



We put 
Then 



H =(L- L)p. 



H(s,t,x,y) = l - g {(Tij (s, x) - <Tij (s, y)) 9 ^ 
■ j = • 

d 



i,j=i ~ l J 



(s,t,x,y) 



dxi 

In the following lemmas k - fold convolution of H is denoted by H^ k \ The following 
results are taken from Konakov and Mammen (2000). 

Lemma 1. LetO<s<t<T. It holds 

oo 

p(s, t, x, y) = ^2p® H (r) (s, t, x, y). 

r=0 

Lemma 2. Let < s < t < T. There are constants C and C\ such that 

\H(s,t,x,y)\ < Cip"Vc, P (y - x) 
\p0H^(s,t,x,y)\ < C^^^y-^y-x), 

where p 2 = t — s, (pc,p( u ) = P~ d 0c(w/p) and 

exp{-C \\u 2 \\) 



<Pc(u) 



J exp(— C \\v 2 \\ dv) 



4 Parametrix method for Markov chains. 

For any < jh < T, x,y G M. d we consider an additional family of "frozen" Markov 
chains defined for jh < ih < T as 

(7) X i+lth = X ith + m (ih, y) h + VhC i+lth , X jth = xeR d , j < i < n, 

9 



where £ n ,h is an innovation sequence such that the conditional density of £,i+i,h 

given X{ t h = Xi,...,X ^ = xq equals to q(ih,y,-). Let us introduce the infinitesimal 
operators corresponding to Markov chains (JTJ) and (JJJ) respectively, 



L h f(jh,kh,x,y) = h 1 IJ p h (jh,(j + l)h,x,z)f((j + l)h,kh,z,y)dz 

-f((j + l)h,kh,z,y)), 
L h f (jh,kh,x,y) = h" 1 (^j ' f? h (jh,(j + l)h,x,z) f ((j + l)h,kh,z,y)dz 

-f((j + l)h,kh,z,y)), 

where p\ (jh,fh, x, •) denotes the conditional density of Xj^h given Xj^ = x. As before 
for any z denote Ph(jh, j'h,x, z) = p\(jh, j'h,x, z) , where the variable z acts here 
twise: as the arument of the density and as defining quantity of the process X^. For 
technical convenience the terms / ((J + 1) h, kh, z, y) on the right hand side of Lhf and 
Lhf appear instead of / (jh, kh, z, y) . 

In analogy with the definition of H we put, for k > j, 

E h (jh, kh, x, y) = (L h - L h J p h (jh, kh, x, y) . 

We also shall use the convolution type binary operation ®h '■ 

fc-i . 

9®hf (jh, kh,x,y) = S^h I g (jh, ih, x, z) f (ih, kh, z, y) dz, 

M jRd 

where < j < k < n. Write g ® h Hf ] = g and g <$ h H { h r) = (g <$ h H^ 1 ^ ® h H h 

for r = 1, ...,n. For the higher order convolutions we use the convention Yl\=j = f° r 
I < j. One can show the following analog of the "parametrix " expansion for ph [see 
Konakov and Mammen (2000)]. 

Lemma 3. Let < jh < kh < T. It holds 

fc-j 

p h (jh,kh,x,y) = ^2p h ®hH { h r) (jh,kh,x,y), 

r=0 

where 

Ph(jh,jh,x,y) =p h (kh, kh,x,y) = 5(y - x) 
and 5 is the Dirac delta symbol. 



5 Bounds on ph — p based on Edgeworth expansions. 

In this subsection we will develop some tools that are helpful for the comparison of 
the expansion of p (see Lemma 1) and the expansion of ph ( see Lemma 3). These 
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expansions are simple expressions in p or ph, respectively. Recall that p is a Gaussian 
density, see (P), and that ph is the density of a sum of independent variables. The 
densities p and ph can be compared by application of the classical Edgeworth expan- 
sions. This is done in Lemma 5. This is the essential step for the comparison of the 
expansions of p and ph- For the proof of Lemma 5 we need one additional lemma 
which is a technical tool for the further considerations. To formulate this lemma we 
need some additional notations. Suppose X e R d be a random vector having a density 
g(x),x &M d , EX = 0,Cov(X, X) = E, where £ be a positively definite d x d matrix . 
Denote A = \\dij\\ = IT 1 / 2 and let Xv(Z) De a cumulant of the order v = (z/i, u d ) of 
a random vector Z e M d , (f>(x) denotes a function in M. d such that D v x ${x) exist and 
continuous for \v\ =4, and A^ 1 = ||a iJ '|| = S 1 / 2 . 

Lemma 4. The following relation holds for s = 3 and for s = 4 



where z = Ax. 

Proof of Lemma 4. For \u\ — 3,u — {v\, v d ), each cumulant \ V {AX) is a 
linear combination of Xfi(X) with = 3 and with coefficients depending only on a^. 
It follows from the following relation 



Xu(AX) = fi u (AX) = / (a u xi + ■■■ + ai d x d y x x ... x (a n xi + ■•• + a ld x d ) Ud q(x)dx. 



Analogously, from the usual differentiation rule of a composite function and from the 
relation (f>(z) = (fi(Ax) , x = A~ x z, it follows that D u z (j)(z) = D u z <p(Ax) is a linear 
combination of D^^Ax) with coefficients depending only on a l K As a result of such 
substitutions we obtain that 





x 





X 




1 d 



fii\...fi d \ 



i 



fll [id / 
J a lj ■■■ a dj a il a kqXfM+ei+e q {^- ) 



' {i^k} l,q=l 1^1 = 1 
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X 



7 y -i-(^Y ; ...(^ d )^a !i 'a fc ''i):' +C, ' +e '>(iX) 

l',q'=l \fj/\=l rl ^ d 



where ( Y2u^j^k}) denotes the sum over all different pairs (triples) of i,j G 

{1, 2, d} ( of k G {1, 2, d}) and e« G M. d denotes the vector whose i— th coordi- 
nate is equal to 1 and other coordinates are zero. Collecting the similar terms in the last 
equation we obtain that for v = 3e k , v' = 3q the coefficient before Xu(X)D^ (j)(AX) 
is equal to ^(a lk a 11 + ... + a dk a ld ) 3 = -^5 kl , for v = e q + 2e r , u' = e/ + 2e n , q ^ r, the 
coefficient before x v {X)D v x <p(AX) is equal to ^(ai q a 11 + ... + ad q a ld ) (a\ r a nl + ... + 
^dr^ nd ) 2 = ^jSqidm, in particular , for I — n the last expression is equal to zero. 
For v = e q + e r + e n , v' = e q > + e r > + e n i q ^ r,q ^ n,r ^ n, the coefficient 
before Xv(X)D^ <f>(AX) is equal to {a lq a q>1 + ... + a dq a q ' d ) x (ai r a r ' 1 + ... + a dr a r ' d )x 
(ai n d n>1 + ... + a dn a n ' d ) = 5 qq <5 rr i5 nn > . This proves lemma for \v\ = 3. The proof for 
\u\ = 4 is quite similar. For this case we use the relation which enabes to express a 
cumulant Xu(AX) as p, v (AX) plus a second order polinomial of the moments fi u >(AX), 
\v'\ = 2. A necessary correction term for n v {X) to get a Xv{X) comes from the deriva- 
tion of D z (j)(z). This completes the proof of the lemma. 

In Lemma 6 bounds will be given for derivatives of ph. The proof of this lemma 
also makes essential use of Edgeworth expansions. The following lemma is a higher 
order extension of the results in Section 3.3 in Konakov and Mammen (2000). Denote 

k-1 k-1 

(8) A*i,fc(j/) = h ^2 m ( ih ' v)> V hk(y) = h y a ( ih i y)- 

i=j i=j 

Lemma 5. The following bound holds with a constant C for v = (y\, ...v p ) T with 
< \v\ < 6 

\D u z p h (jh, kh, x, y) - D u z p(jh, kh, x, y) 
—\Z~hD 1/ z Tii(jh, kh, x, y) — hD z n 2 (jh, kh, x, y) 



<ch*" P -*c;-M(v- 



X) 



for all j < k,x and y. Here D z denotes the partial differential operator of order v 

—ill 

with respect to z = V- k {y){y — x — fij <k (y)). The quantity p denotes again the term 
p— [h(k — j)] 1 / 2 and the functions 7?i and 7r 2 are defined in (J3J and (JEJ). We 
write Cp(') = P~ P C h (-/p) where 



f[l+ Wz'fj^dz' 



Proof of Lemma 5. We note first that Ph(jh, kh, x, •) is the density of the vector 

k-1 

i=j 
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where, as above in the definition of the "frozen" Markov chain Y n , £i+i ( h is a sequence 
of independent variables with densities q(ih,y,-), Pj, k (y) = ^2t=j hm(ih, y) . Let fh(-) 
be the density of the normalized sum 

fe— 1 

h 1/2 [Vj, k (y)r 1/2 J2l + i, h . 

i=j 

Clearly, we have 

p h (jh, kh,x, •) = det [Vj, k (y)r 1/2 fn{{v J)k (y)Y 1 ' 2 [• ~ x - mMW- 

We now argue that an Edgeworth expansion holds for fh- This implies the following 
expansion for Ph{jh, kh, x, •) 

(9) p h (jh,kh,x,-) 

S—3 

= det [v hk (y)r l/2 [J> - j)- r/2 P r (-<t> : {xeA)WMr 1/2 [• - * - hAv)]} 

r=0 

+[k - j]- {s - 2)/2 0([l + \\{[V j>k (y)r 1/2 [- - * - Hk(y)]}\\ S r 1 )] 

with standard notations, see Bhattacharya and Rao (1976), p. 53. In particular, P r 
denotes a product of a standard normal density with a polynomial that has coefficients 
depending only on cumulants of order < r + 2. Expansion (JHJ) follows from Theorem 
19.3 in Bhattacharya and Rao (1976). This can be seen as in the proof of Lemma 3.7 
in Konakov and Mammen (2000a). 

It follows from (JH]) and Condition (A3) that 

\ph{jh, kh,x,y) -p(jh, kh,x,y) -h 1/2 9i(jh, kh,x,y) - hTT 2 (jh, kh,x,y)\ 



(10) 
where 



<C0 2 p-^%-x), 



p(jh, kh, x, y) = det [V jjk (y)\ (2tt 



1/2 , 9(jr N-p/2 



cx p{-2^ ~ x ~ vjM) t MM] (y~ x - VjM)}, 



9i(jh, kh, x, y) = -p 1 det [V jtk (y)] 1/2 ^2 



-i/ 2 \- ^^ D ^\ [VjM r 1/2 iv-x- HM) 



["1=3 



7T2(jh, kh, x,y) = p 2 det [Vj 



j,k 



-1/2 



E 

M=4 



Xuj, k (y) 



D u jm k ^- 1/2 



(y - x - n jtk (y))j 



+ E 



x»,j, k (y) 



2 |" I/! 

I"l=3 
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where X v ,j,kiv) = 7TJ Yh=} Xvj,k,i(v)> X»,j,k,i(y) = f-th cumulant of p [V jtk (y)] 1/2 
p\ v \ x {z/— th cumulant of [Vj^iy)]' 1 ^ 2 £i+i,h}, and D u z (j)(z) denotes the v— th derivative 
of 4> with respect to z = \Yj,k{y)]~ 1 ^ 2 (y — % — ^j,k(y)) • It follows from the (conditional) 
independence of = j, k - 1, that X v ,j,k(y) = k~j h ~ M/2 x Xu(AX), where 

^4 = h^[V hk {y)}~ 1/2 = E-V^E = CW(X,X),X = Et/W By Lemma 4 for 



3,4 



|i/|=s kl=s 

\v\=s 

(ii) = (-i) V £ ^^ DV Mv h k{y)V 112 (y - * - m*(v))), 



Z/! 



where we put 0h(z) = 0(/i~ 1 ^ 2 z), ^(-X") = Ei=j Xu(^h,y). It follows from (JTTJ) and 
the condition Bl that up to the error term in the right hand side of (fTUl) the functions 
7Ti and 7?2 coincide with the functions 7F1 and tt 2 given at the beginning of Section 4. 
For v = the statement of the lemma immediately follows from (fT0|) . For v > one 
proceeds similarly. See the remark at the end of the proof of Lemma 3.7 in Konakov 
and Mammen (2000). 

The following lemma was proved in Konakov and Molchanov (1984) (Lemma 4 on 
page 68). 

Lemma 6. Let L(d) be the set of symmetric matrices, and let D\+ \- C L(d),0 < 
A~ < A + < 00, be the open subset distinguished by the inequalities A G D\+ \- <^> 
X~I < A < X + I; A = A(A) is a solution of the equation A 2 = A with A = A T and 
ajj(A) the corresponding matrix elements. Then for any k,l < d,i,j < d and 
A C D\+ \- we have that 



daij(A) 



dX 



kl 



1 

< 



2VA 



From Lemma 5 we get the following corollary. The statement of the next lemma 
is an extension of Lemma 3.7 in Mammen and Konakov (2000) where the result has 
been shown for < |fe| < 2, a = 0. 

Lemma 7. The following bounds hold: 

\D;D b x p h (jh,kh,x,y)\ < Cp-\ a ^ h k s p ~\ a \y 



x 



for all j < k, for all x and y and for all a,b with < \a\ + \b\ < 6. Here, 
p = [(k — j)h] l l 2 . The exponent S has been defined in Assumption A3. 

Proof of Lemma 7. If A = \\aij\\ and B = and elements a,ij(B) are 

smooth functions of b^i then an inequality |§4| < C will mean that 



r)b,. 



< c 
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for all 1 < i,j < d,l < k,l < d. To obtain the assertion of the lemma we have 
to estimate the derivatives DyD^z, where z = VjJ (y)(y — x — ^j,k(y))- Note that 
z = z(V~ k 1/2 , /J, j: k , x, y), where = V^ k 1/2 {y) and fi jtk = fij,k(y)- It follows from the 



conditions (Bl) and 
(12) 



that 

d V>3,k{y) 



dy 



It follows from Lemma 6 that 



(13) 



<Cp 2 



dV j>k (y) 



1 k 



dy 



< C. 



<Cp 2 . 



From inequalities (3.16) in Konakov and Mammen (2000) and from the representation 
of an inverse matrix in terms of cofactors divided by determinant we obtain that 



(14) 



< c P - 



(15) 



From (|r^jl -<[n jl and from the differentiation rule of a composite function we get 

9Vr k 1/2 (y) 



dy 



< C. 



Inequalities (3.16) in Konakov and Mammen (2000), (fTBTl and the differentiation 
rule of a composite function imply 



(16) 



dz 



dpij 



3,k 



< C -, 

P 



dz 



{y-x- iijjk(y)), 



dz 



dy 



< c -. 



Inequalities (3.16) in Konakov and Mammen (2000) also imply 

dz 



(17) 



dx 



C 
< —. 

P 



The assertion of Lemma 7 for a = e», b = ej, 1 < i, j < d, follows from Lemma 5 and 
from ((HJ), (fTTj) . For other values of a and b one has to repeat this arguments. 



6 Bounds on operator kernels used in the parametrix 
expansions. 

In this section we will present bounds for operator kernels appearing in the expansions 
based on the parametrix method. In Lemma 8 we compare the infinitesimal operators 
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Lh and Lh with the differential operators L and L. We give an approximation for the 
error if, in the definition of Hh = (Lh — LhjPh, the terms Lh and Lh are replaced by L 
or L, respectively. We show that this term can be approximated by K h + M h , where 
Kh — (L — L)ph and where Mh is defined in Remark 1 after Lemma 8 . Bounds on 
Hh, Kh, and Mh are given in Lemma 9. These bounds will be used in the proof of our 

~ (r) 

theorem to show that in the expansion of ph the terms ph ®h H h can be replaced by 
p h Oh (K h + M' h Y r \ M' h is defined in Lemma 9. 

Lemma 8. The following bound holds with a constant C 

\H h (jh, kh,x,y) - K' h (jh, kh, x,y) - M' h (jh, kh,x,y) - R h (jh, kh, x,y)\ 
< C0 2 p-\ s p (y - x) 

with (p as in Lemma 5 for all j < k, x and y . For j < k — 1 we define 

K' h (jh, kh, x, y) = (L- L)X(x), M' h (jh, kh, x, y) 

= M hA (jh, kh,x,y) + M h , 2 (jh, kh,x,y) + M' h3 (jh, kh,x,y) 

M K1 (jh,kh,x,y) = h 1 ' 2 Dx l{ iXu(jh,x) - Xu(]h,y)), 



M=3 



M h , 2 (jh, kh, x,y) = h^2 — — (Xu(jh, x) - Xu(jh, y)), 



M=4 

M^(jh,kh,x,y) = ^(Ll-L 2 )\(x) 



D v \(x) 

R h ( ] h,kKx,y) = 2 Y,^T 1 Y. 

u r [m r (jh,x)p, u - er (jh,x) - m r (jh,y)fi u _ er (jh,y)] 

H=4 ' r=l 

+ 5 ^^Y^( m k(jKx)-rn k (jh,y)) \v k j q(jh, x, 6)1^(6) 

\v\=5 V ' k=l ^ J 



X 



A- 

Jo 



uYD v \(x + uh(6))du 



d9 



+ / q(jh,x,9)h v (9) 



1 — u) uD v+ek \(x + uh(9))du 



d9 



+h 2 ^2 ~^T~ v } - N (v,v')[™ u '(jh,x)n l/ - l// (jh,x) - m u ' (jh,y)p, v _ v ,(jh,y)\. 



| z,|=4 



M=2 



Here L* is defined analogously to L but with the coefficients "frozen" at the point x, 
e r denotes a p - dimensional vector with r — th element equal to 1 and with all other 
elements equal to 0. Furthermore, for |z/| = 4, \u'\ =2 we define 

N(v, v) = 2xK !=1 ]+x[(^-^)!=i]-2 

where x(') means an indicator function. We put m(x) v = mi(x) Vl ■ ... ■ m, p (x) Vv and 
m(x) u — 0, v\ — and p. v (x) = if at least one of the coordinates of v is negative. 
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We define also the following functions 

K x ) = Ph(U + l)h,kh,x,y), 
h{6) = m(jh,y)h + 9h 1 l 2 . 

1 /2 

Here again p denotes the term p = [h(k — j)] . For j — k — 1 and I = 1, ...,3 we 
define 

K' h (jh,kh,x,y) = M h>1 (jh,kh,x,y) = M h>2 (jh,kh,x,y) = M' h3 (jh, kh,x,y) = 0. 

Proof of Lemma 8. As in the proof of Lemma 3.9 in Konakov and Mammen (2000) 
we have 

H h (jh, kh, x, y) = H l h (jh, kh, x, y) - H 2 h (jh, kh, x, y), 

where 



(18) Hl(jh,kh,x,y) = h- 1 j q(jh,x,9)[X(x + h(9)) - X(x)]d9 

(19) Hl(jh,kh,x,y) = h' 1 J q(jh,x,9)[X(x + h(9)) - X{x)]d9, 

h{9) = m(jh,x)h + 9h 1/2 , 
h{9) = m{jh,y)h + 9h 1/2 . 

For [X(x + h(9)) - X(x)} and [X(x + h{9)) - X{x)\ in (IT8l). (fTTil) we use now the Tailor 
expansion up to the order 5 with the remaining term in integral form. To pass from 
moments to cumulants we use the well known relations (see e.g. relation (6.11) on 
page 46 in Bhattacharya and Rao (1986)). After long but simple calculations we come 
to the conclusion of the lemma. 

Remark 1. We show now that the function K' h {jh,kh,x,y) + M' h3 (jh,kh,x,y) 
in Lemma 8 is equal to K h (jh, kh, x,y) + \{L 2 — 2LL + L 2 )X{x) + Mjl 3 (jh,kh,x,y) 
where 

M'^jh, kh, x, y) = -h 2 m " U ^ y) (L - L)D»\(x) 

(20) -3V f (l-5) 2 d5 [ q(jh,y,9)^^(L-L)D^X(x + 5h(9))d9. 

Thus in Lemma 8 we can replace K' h (jh, kh, x, y)+M' h (jh, kh, x, y) by K h {jh, kh, x, y)+ 
M h (jh,kh 2 x,y) where K h (jh,kh,x,y) = (L - L)p h (jh,kh,x,y), M h (jh,kh,x,y) = 
^{Ll-2LL+L 2 )X{x)+Ml M'l = M h)l (jh, kh, x, y)+M h , 2 (jh, kh, x, y)+M'l 3 {jh, kh, x, y) 
and 

max{\M' h (jh, kh, x,y)\, \M h (jh, kh, x,y)\} < Cp^Cpiy - x), 
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p 2 = kh — jh. To show this we note that 



PhUh,kh,x,y) = J q(jh,y,0)\(x + h(0))ffl 
where h(9) = m(jh,y)h + h 1 l 2 9. From the Tailor expansion we get 

p h (jh,kh,x,y) = X(x) + hLX(x) + h 2 J2 m „, fl^r) 



M=2 

+3 J2 I (l-$) 2 d6 [ q(]h,y,9)^I^D^\(x + 5h(9))d9 
M=3 Jo J /i - 

and, hence, 

K' h (jh,kh,x,y) = K h (jh,kh,x,y) + (L - L)[X(x) - p h (jh, kh, x,y)} 



//! 



(21) 

Note that 



K h (jh, kh, x, y) + h(L 2 - LL)X(x) + M% 3 (jh, kh, x, y). 



h 



h{L 2 - LL)X{x) + M' h3 {jh,kh,x,y) = h[L 2 - LL)X{x) + -[L\ - L 2 )X[x) 



h 



(Ll - ILL + L 2 )X{x) 



and from the definitions of the operators L, L and L* and Lipschitz conditions on the 
coefficients m{t,x) and a{t,x) we obtain that 



(22) 

Analogously, we have 



h 



(Ll-2LL + L 2 )X(x) 



<Chp- 3 ( p (y-x). 



(23) 
(24) 



h ^ m ^y\ L - LWX {x) 



H=2 



//! 



<Ch 2 p - 3 c P (y-x), 
h{9Y 



3 V [ {\-5fd5 [ q(jh,y,9)^^(L- L)D»X{x + 5h{9))d9 



|A*I=3 

<Ch 3 / 2 p-%(y-x). 



Now (|2T|l - (f24|) imply the assertion of this remark. 
Lemma 9. The following bounds holds: 

n n 

Ph ®h {K h + M h + i4) (r) (0, T, x,y)-J2 Ph ®h {K h + M h )^ (0, T, x, y) 
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r=0 



r=0 



(25) 



<C(E)hn- l l 2+ %: Uy-x), 



where lim e >oC(e) = +00. 

Proof of Lemma 9. For r = 1 we will show that for any e > with p 2 = kh 

(26) \p h ® h (K h + M h + R h )(0, kh, x, y) - p h ® h (K h + M h )(0, kh, x, y) \ 

= \p h ® h R h (0,kh,x,y)\ < Ch^ 2 - e (kh)-^ 2+e B(^,e)C^y - x). 

Clearly, for an estimate of ph <S)h Rh(0, T, x, y) it suffices to estimate 

fc-2 . 

I, = tf' 2 J^h Ph(0, kh, x, z)(f(jh, z) - f(jh, y))D v z p h ((j + l)h, kh, z, y)dz 

j=0 J 
for v with \v\ =4, and 



fc-2 



I 2 = h 2 J2h[p h (0,jh,x,z)(f(jh,z)-fUh,y))fq(jh,z,9)h^(e) 

x f (l-ufD v z \(z + uh{d))dud6dz 

for v with \v\ — 5, 1 < k < d. Here f(t,x) is a function whose first and second 
derivatives with respect to x are continuous and bounded, uniformly in t and x. After 
integration by parts we obtain for 1 < I, s < d 



-0 2 J^h D e z «p h (0, jh, x, z)(f(jh, z) - f(jh, y))D:- e "p h ((j + l)h, kh, z, y)dz 
3=0 J 

fc-2 „ 

+/1 3 / 2 J^h D e /p h (0, jh, x, z)Dt"f(jh, z)D:~ e ^p h ((j + l)h, kh, z, y)dz 
3=0 J 

fc-2 „ 

+0 2 Y,h p h (0, Jh, x, z)D e z ^f(jh, z)D^ e ^p h ((j + l)h, kh, z, y)dz. 
1=0 J 



3=0 
Hence, 

fc-2 



(27) \h\ < Ch*l*Y, h mirh - h M v-*) < Ch 3 / 2 ^(kh)-^B(-,e)C%(y-x). 

VJh{kh-jh) 2 V1 



j ..... ., . , ,, 1 

In the same way, we obtain after integration by parts with 1 < Z, 5 < d 

k—2 1 

I 2 = -h 2 J2h[ (l-u) 4 du [ d6(m(jh,y)h 1/2 + 6y~ ei 



D?p h (0, jh, x, z){f{jh, z) - f(jh, y)) 
xq(jh, z, 9)D^ e 'p h ({j + l)h, kh, z + uh{6),y)dz 
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^ 2 pi p 

+h 2 J2 h / {l-ufdu / d6(m(jh,y)h 1/2 + ey- £l 
x f D e /[p h (0,jh,x,z)D^f(jh,z)q(jh,z,e)} 



xD:- e >~ e °p h ((j + l)h, kh, z + uh(6),y)dz 

J2 h {^-ufdu / de{m{jh, y )h}i 2 + ey~ ei 

xp h (0,jh,x,z)(f(jh,z) - f(jh,y)) 

xD e z l q(jh, z, 9)D u z - £l p h ((j + l)h, kh, z + uh(6),y)dz. 



It follows from this equation that 

(28) |/ 2 | < Ch^ikhr^B^eK^y-x). 

Claim (|2*T|) follows now from (|2*T|) and (EHJ) - For r > 2 we use the identity 

(29) p h ® h (K h + M h + i4) (r) (0, T, x, y) - p h ® h (K h + M h ) (r) (0, T, x, y) 

= [p h ® h (K h + M h + i4) (r_1) - Ph ® h (K h + M h f-^] ® h (K h + M h )(0, T, x, y) 

+p h ® h (K h + M h + R h f^ ® h R h (0, T, x, y) 
= 1 + 11. 

For r = 2 we obtain from (|2Tjl and from the simple estimate | (Kh + Mh) (jh, kh, z,y)\ < 
Cp 2 1 Cp 2 (y — z) with p\ = kh — jh, that 

|/| = \\p h ^ h (K h + M h + R h )-p h ^ h (K h + M h ))^ h {K h + M h )(0,kh,x,y)\ 

< C 2 h^B{\,e)Y,h{ 3 h)- l ^ {kh _ jh) -i/2 f c ^ {z _ x)( -s^ y _ z)dz 

1 3=0 J 

< C 2 2 ~^B{ 1 -, e)B{ 1 -, e + l -){kh)% s p {y -x),p 2 = kh. 

For r > 3 we obtain by induction 

|/| = \[p h ® h (K h + M h + R h Y r -V -p h ® h ( Kh + M h Y r -V] 
® h (K h + M h )(0,kh,x,y)\ 

< C h^B{\,e)B{\, e + \)...B{\, e + Lll) W ^% 5 (y " x) 

(30) < V{e)h^ [ °^1 2 }J (khy^-^iy ~x),P 2 = kh. 
To estimate II we use the following estimates 

(31) \D a y D%(jh,kh,x,y)\ <Cp-^ b ^(y-x),D%(jh,kh,x,x + v)<C^(v), 

(32) \D b x {K h + M h + R h )(jh, kh, x + v, x)\ < Cp-'C^v). 
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The inequalities (|31f and (|32j) are obained by using the same arguments as is the proof 
of Lemma 7. Using these inequalities and mimicking the proof of Theorem 2.3 in 
Konakov and Mammen (2002) we obtain the following bounds for r > 

\D b x D°p h ® h (K h + M h + R h ) {r) (0, kh,x,y)\ 

< C7(A*)-W-W^B(i \)B{1, - x) 

(33) < M(tt)- H * r <- |a| (!/-4 

Inequality |3T?| allows us to estimate II = [ph®h(Kh + Mjl + Rh)( r ~^]®hRh(0, kh, x,y). 
For this aim it suffices to estimate 

h 3/2 J2 h J (K h + M h + R h f-V](0,jh,x,z) 

(34) xD»p h ((j + l)h, kh, z, y)(f(jh, z) - f(jh, y))dz 
for r > 2 , \v\ = 4, and 

n— 2 „ 

h / [p h ® h (K h + M h + R h f-V](0, jh, x, z)(f(jh, z) - f(jh, y)) 

3=0 J 

(35) x [ q(jh,z,9)h u - e '(9) [ (1 - u) 4 D u p h ((j + l)h, kh, z + uh{9),y)dud9dz 



for r > 2, \v\ — 5, 1 < I < d. Here f(t,x) is a function whose first and second 
derivatives with respect to x are continuous and bounded, uniformly in t and x. The 
upper bound for (p?Ijl follows from (j3TS|) by integration by parts, exactly in the same 
way as it was done to obtain the upper bound for Ii (see ([27]) ). This gives an estimate 



ti" 2 



J2 h J &h®h (K h + M h + R h )^](0,jh,x,z) 
xD»p h ((j + l)h, kh, z, y)(f(jh, z) - f(jh, y))dz\ 



(36) < r( g )^ 2 - ig^C (^)^- 2 )/ 2 c P 5 (y - x). 

The upper bound for (|35l) also follows from (|33| by integration by parts in the same 
way as it was done to obtain an upper bound for I 2 (see (|2lfl) ). This gives for (j3~51) the 
same estimate as in (j3~El) and, hence, 

(37) |//| < CT{e)h^ [C l^ T {khf^y\ s p {y - x). 
The assertion of the lemma follows now from (|2*7j) . (JSHJ), (|SIH) and (pTTjl . 

Lemma 10. Ze£ A(s,t,x,y), B(s,t,x,y),C(s,t,x,y) some functions with the 
absolute value less than C(t — s) _1 ' 2 C /^(y — x) for a constant C. Then 

oo oo 

A ®h (B + C) (r) (ih, jh, x,y)-J2 A ®h B {r) (ih, jh, x, y) 

r=0 r=0 

oo 

= l A ®h $] ®fc [C ® h $] (r) (i/i, jh, x, y) , 



r=l 
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where $ = E^o fiW - 

Proof of Lemma 10. Under the conditions of the lemma all series are absolutely 
convergent. The assertion of the lemma follows from linearity of the operation ®h and 
after permutation of the summands in the absolutely convergent series. 



7 Proof of the main result. 

We now come to the proof of Theorem 1. The main tools for the proof have been 
given in Sections 4 and 5. From Lemmas 1 and 2 we get that 

n 

p(0, T, x, y) = ® # (r) (0, T, x, y) + o{h 2 T)^ c ^{y - x). 

r=0 

With Lemma 3 this gives 

(38) p(0, T, x, y) - Ph (0, T, x, y) = T x + ... + T 7 + o(h 2 T)<P c ^ (y - x), 

where 

n n 

Tl = P ® H W (0, T, x, y) - P ®h H (r) (0, T, x, y), 

r=0 r=0 

n n 

T 2 = J2P®h H^ r \0, T,x,y)-J2P ®h (H + M£ + Vh^ ) M (0, T, x, y) , 

r=0 r=0 

n 

T 3 = ^2p® h {H + Mjl + VhNjW ( 0) T , x, y) 

r=0 

n 

- P ®h (H + M h + VhN ± ) W (0, T, x, y) , 

r=0 

n n 

T, = Y,P®h{H + M h + VhN,) W (0,r,i,3/)-^p ®, + M fc ) (r) (0, T, x, y) , 

r=0 r=0 

n ra 

T 5 = P ®h (K h + M h Y r \0,T,x,y)-Y, Ph ®h (K h + M h ) W (0, T, x, y) , 

r=0 r=0 

n n 

T 6 = ]T p h ® h (K h + M fc ) « (0, T, x,y)-£ ®, (X, + M h + R h ) < r > (0, T, x, y) , 

r=0 r=0 

n n 

T 7 = ® h (tf ft + M h + W(0, T, x, y)-J2Ph ®h H ( h r) {0, T, x, y). 

r=0 r=0 

Here we put Ni(s,t,x,y) = (L — L)tti (s, t , x, y) . 

We now discuss the asymptotic behavior of the terms Ti, ...,T 7 . 
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Asymptotic treatment of the term T\. Using Remark 1 in Konakov (2006) we get that 

T x = ^[p® h (L 2 -2LL + L 2 )p® h $)](0,T,x,y) 
h, 



(39) 

where for any < e < 1/2 



+ -[p ® h (L' - L')p ® h $](0, T, x, y) + R T (0, T, x, y), 



(40) 



\R T (0,T,x,y)\<C(e)(hn 



-1/2+6 



hVT)<t) Cy/T (y 



x 



&(s,t,x,y) = Y^^=o H { s i t-i x i v) ■ Here, the summand H^\s,t,x,y) was introduced 
to simplify notations. We define g ®h H®>(s,t,x,y) = g(s,t,x,y) for a function g. 
Note, that in the homogenous case it holds that (Tij(s, x) = aij(x),rrii(s,x) = rrii(x) 
and that the second summand in (J3*9"J) is equal to 0. 

Asymptotic treatment of the term T 2 . We will show that with a positive constant 
5 > 



(41) 



T 2 -3j2P®hH (r \0,T,x,y) + J2P®h(H + M K1 + VhN^ (0, T, x, y) 



r=0 



r=0 



+ J2p®h(H + M h ^ (0,T,x,y) + J2P®h(H + Mj 3 ) W (0, T, x, y) 



r=0 



r=0 



<Chn~ 5 CsT{y-x). 



It suffices to discuss the case r > 2 because for r = 1, 2 the left hand side of (|4T|) is 
equal to zero. For r > 2, (pTl) immediately follows from the following bounds 



(42) 



(43) 



p ® h (H + Ml + VhN t ) W(0, T, x, y) 

-p ®, (i? + M M + M h)2 + VWVi) (r) (0, T, x, y) 
-\p ® h (H + Ml,) « - p ® fc fT«](0, T, x, y) | 

p® h (H + M h>1 + M K2 + v^iVi) (r) (0, T, x, y) 

-p ®, (if + M M + V^JVx) W(0, T, x, y) 
-\p ® h (H + M h>2 ) W - p ®, ffW] (0, T, x, y) | 

for all sufficiently small e > with a function C that fulfils lim £ , C(£) = +oo. We 

will prove first the bound (|4*2l) . Denote the expression under the sign of the absolute 
value in (|4~2|) by T r . Note that r = Ti = 0. For r > 2 we make use of the following 
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recurrence formula 

(44) r r = r r _! ® h h + 



(H + M'l + VhNiY r ~ 



i) 



(H + M M + M K2 + VhNx) 



(r-l) 



{Ml + y/hNj 



+ \p ® h (H + M h>1 + M h>2 + VhN^-V -p® h (H + MZ >3 )<r-V 
= 1 + 11 + III. 

We start with the estimation of II. First we will estimate 



h M. 



h,3 



(45) 



{H + Ml + VhNJ 



(r-l) 



P®h 



{H + M h>1 + M K2 + y/hN 1 ) 



(r-l) 



For r = 2 we have p ®h M^ 3 (0, kh, x, y). It follows from (|20|) that it is enough to 
estimate 

fc-2 „ 

(46) J 1 = h 2 J2h p(0,ih,x,v)(f(ih,v) - f{ih,y))D v v p h ((i + l)h,kh,v,y)dv 



for \v\ = 4 and 



fe-2 



(47) J 2 = /rV 2 ]T/i f p(0,ih,x,v)(f(ih,v)-f(ih,y)) f q(ih,- 
xD^ +ei+e "p h ((i + l)h, kh, v + 5h(6),y)d5d6dv 



(1-8)' 



for \u\ = 3. Here f(t,x) is a function for which D^f(t,x) is bounded uniformly in (t,x) 
for \u\ = 0, 1,2,3. An estimate forJi follows from (J2"?jl . This gives 

(48) \J X \ < Ch^ikhr^B^eK^y-x). 

The estimate for J 2 can be obtained analogously to the estimate of I 2 (see (|2*Hll ). 
Integrating by parts we get 

J 2 = h 3/2 J2 h / (l-5) 2 ^ d6-6 u D e v l+e "[p(0,ih,x,v) 

o_n ^0 J J 



3=0 



x (/(zTi, u) - /(z/i, y))g(z/i, u, 9)}D v v p h {{i + l)h, kh, v + 8h(9),y)dv. 



The derivative 



ih, x, v)(f(ih, v) - f(ih, y))q{ih, v, 6)] 



is a sum of 9 summands. Integrating by parts once more for summands which contain 
Dffi(0,ih,x,v) with < 2, we obtain 



(49) 



|J 2 | < Ch*^*^{y-x)Jm\\n*(hV- 



3)/2 || fl ||S-3 



+ l)d6 



3=0 



(ihy-z (kh-ihy- 



< Ch^B^e^khf^C^iy-x) 
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for any e G (0, 1/4). It follows from P8|) and (JlHj) that for r = 2 ijlHJ) does not exceed 

r > 3 we use the recurrence relation 

(50) p® h (H + Ml + VhNi)^'" 1 ^ -p® h {H + M M + M hj2 + VhN^ 1 ^ 

= p® h (H + Ml + VhN^-V -p® h {H + M M + il4, 2 + VhNif-V 

® h (H + MI + VhN,) + [p ® h (H + M M + M M + V^V^" 2 )] ® h M£ 3 
= J' + II'. 

From ([50j) we obtain for r = 3 

fc-2 

|/'| < Ch 3 / 2 -^B(e,e)C s ^ h (y-x)Y,h^hy- 1 (kh-zh)- 1 / 2 



i=0 



< Ch^B(e,e)B(l,s)(khry 2 C S ^ h (y 



x) 



For the estimate of IV we use the following estimates 



(51) 
(52) 



D a v D b x (H + M M + M M + VhN x )(jh, kh, 



< Cp 



-l-|a|-|6| 



C P (v-x), 



D b (H + M M + M K2 + VhN{){]h, kh, x,x + v) < Cp' 1 ^ - x). 



To prove (|5Tf it is enough to get the corresponding estimates for summands in M h;1 , M h>2 
and \/hNi 



(53) 
(54) 
(55) 



h^ 2 D a v D b x [D v x p h {{j + l)h,kh,x,v){f{jh,x) - f{jh,v))\ 

< Cp-^-^Cpiv - x) for \v\ = 3, 

hD a v D b [D»p h ((j + kh,x,v)(f(jh,x) - f(jh,v))] 

< Cp-^-^Cpiv - x) for \v\ = 4, 

[Z^ +e ^((j + x, v)p 2 (/(jfc, x) - /(jfc, v))] 

< Cp-^-^Cpiv - x) for \u\ = 3, 



for a function f(t, x) with \a\ + \b\ derivatives w.r.t. x that are uniformly bounded w.r.t. 
t. These estimates are direct consequences of Lemma 7. To prove (|52|) it is enough to 
show the corresponding estimates for summands in M h l , M hj2 and yhNi 

h^ 2 D b [D» x p h {{j + l)h, kh, x, y) \ y=x+v (f(jh, x) - f(jh, x + v))} 



(56) 
(57) 
(58) 



< Cp C p (v — x) for \v\ = 3, 

hD b x [D»p h ((j + l)h, kh, x, y) \ y=x+v (f(jh, x) - f(jh, x + v))] 

< Cp~ l C P {v - x) for \u\ = 4, 

h x l 2 D b x [D^ e ^p h ((j + l)h, kh, x, v) \ y=x+v p 2 (f(jh, x) - f(jh, x + v 

< Cp~ 1 ( p (v — x) for \u\ = 3. 



These estimates also follow from the estimates obtained in the proof of Lemma 7. If 

z (V 3 \k (y)^jAy)^ x ^ y) = v j\k /2 (y)(y - % - mM) then 



(59) 



dz 
dy 



C 

< ~, 
P 



dz 
dx 



C 
< — 

P 
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and for z(V- ^ 2 (x + v),p,j^{x + v), x, x + v) = V- ^ 2 {x + v)(v — fij,k(x + v)) we have 



(60) 



dz 



dx 



< 



dz 



dx 



dz 



dp 



3;k 



dx 



<C(\\v\\ + l). 



Now with the inequalities . we can proceed like in the proof of Theorem 2.3 in 
Konakov and Mammen (2002). This gives the following estimate for r > 3 



(61) 



D a v D b x \p® h (H + M M + M hj2 + y/hN^-^ijh, kh, x, v) 
< C r B(l, I) x ... x B( r -^, l)p r - 2 -^%(v - x). 



Now we denote pi jT = p ®h (H + M^i + Mh,2 + VhNi)^ and po = p. For an estimate 
of Pi >r -2 ®h it suffices to make the same calculations using integration by parts as 
it was done above for J\ and J 2- This gives 

(62) \II'\ < |pi ir _2® fc Mj 3 (0,A;/i,x,y)| 

< Ch**-*B{1, \)B{1 + 1,1) 

x... x 5(1 + \)B( r -^, e)(khy^(^ h (v - x) 



and by induction for r > 3 
(63) 



'\<C r h^B{e,\)B{e + \, \) 



x... x B(e + ~)B(e, £)(kh) s+ ^(^(v - x). 



Comparing (jfi2~|) and we obtain that for r > 3 



(64) 



p ® h (H + Ml + VhN!)^ -p® h 

<C'h^B(e,\)B{e + \,\) 
r-3 1 



(H + Mw + Mw + VhN!) 



(r-l) 



x... x S(e + 



^M)(^) £+ ^Cy^ 



■u — x . 



From (JBIJl we get the following estimate for II 
(65) |//|<C^ 3 / 2 - 2£ 5(e )£ )5(^ 

To estimate III note that the same inequalities (JHH), . (|fiT|) hold for H + M^ 3 , 
namely, 



(66) 
(67) 
(68) 



|L>^(tf + M'l 3 )(jh, kh,x,v)\ < Cp-'-^Uv - x) 
\D b x {H + M£ )3 )(jh, kh, x,x + v)\ < Cp-\ p {v - x) 
\D a v D b x [p® h (H + Ml 3 )^](jh, kh, x, v)\ 

< C r B(l, 1) x ... x b£±±, l)p r - laH %(v - x). 
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To prove (|66|) - (j68|) it is enough to get the corresponding estimates for summands in 
M^ 3 (see (|2fll) ). These estimates can be proved by the same arguments as used in the 
proof of (|5TJ) . (|52~| and (jfiTJl . To estimate J// we have now to estimate p^ r and 
p 2 ,r ®h M'^ 3 where 

Using integration by parts and inequality (j68|l . we obtain for P2,r-i <8>/i the same 
estimate as for pi )T . Cgi^ M^ 3 . It holds for z = 1, 2 

|p 2 ,r ®fe M'h 3{0,kh,x,y) I 

< CT 3 / 2 - 2£ 5(l,i) x ... x B(^,±)B{le)(kh)*TC^y-x). 

Hence it holds for r > 2 



(69) 



\IH\ < \pi t r-i®hM£ 3 (0,kh,x,y)\ + \p2,r-i®hMZ 3 (0,kh,x,y)\ 



< C r h^B{l,\) x ... x B( r -,l)B(^,e)(khy 



From (HU), jnnj) and fIE|) we get for r > 2 

|r r (0, kh, x, y)\ < C r h^B{e, e)B(e, i) x ... x B(e + \){khf +z ^ C m {v - x). 



In particular 



(70) 



\T r (0,T,x,y)\ <h^ 2 - 



2c 



C 



T(2e)T(e+ r -^ 



-T £ 



V — X) 



for any e G (0,1/4) and r > 2. Now we shall estimate the left hand side of (|4*3|) . 
Denote the expression under the sign of the absolute value in (|4~H|) by F r - Note that 
Fo = Fx = 0. For r > 2 we make use of the following recurrence formula 



r, 



r, 



r-l 



® h (H + M M + M M + VhN^-V 



-p® h (H + M M + VhN^-V] ® h (M M + M K2 + VhN, 



+ [p® h {H + M M + VhN^-V -p® h (H + M htl )<- r -Q 

i + 11 + in. 



We start again from the estimation of 

A r -i =p® h {H + M M + M h>2 + VhN^-V -p® h (H + M M + VkN^l 
For r = 2 we have A x = {j>®h M h2 )(0, kh, x, y). It is enough to estimate 

J 3 = h^h / p(0, ih, x, v)(f(ih, v) - f(ih, y))D u u p h ((i + l)h, kh, v, y)dv 
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for \u\ = 4. Analogously to (|27j) we obtain that 



\Jz\ < Ch^kh)-^B(\,eK: u { !f - ., ) 



\A t \ < Ch^ikhy^B^eX^y-x). 



and, hence, 
(71) 

For r > 3 we use the recurrence relation 

(72) A r _i = A r _ 2 ®h {H + Mh,x + Mh,2 + vTtiVi) 

+ \p® h {H + M M + v^iVi) {r - 2) 
= /' + //'. 

From (I7TJ) and (f72~| we obtain for r = 3 



h M h o 



1 fc-2 

|/'| < Ctf-'Bi-rfCfafa -x)J2Hih) £ - 1/2 (kh-ih) 



8=0 



< ^"^(i e)S(i £ + i)(^) £ c5o:(y - x). 



(73) 



'2' 7 v 2' 2 
To estimate //' we use the inequality for r > 3 

D a v D b x [p® h (H + M M + VhN^-^jh, kh, x, v) 



-1/2 



(74) 



This inequality is a direct consequence of (j6Tl) . We have 
(75) \II'\ < Ch^B^e^khYC^-Ay-x). 



Comparing (jZBJ) and (ZED we obtain that \A 2 \ < C 2 h 1 - £ B(±,e)B(±,e + l)( y kh) £ (^ J - 
x). By induction we easily get that for r > 2 

(76) \A T - X {Q,kh,x i y)\<Crh 1 -°B{± 1 e) x ... x 5(1 e + r ^)(khf +I ^ ^{v ~ 
For a bound of A r _i ®h (M^i + M^2 + VhNi) it is enough to estimate 



fc-2 



(77) 



(78) 



(79) 



j 4 = h V*J2h [ A r ^(0,ih,x,v)(f(ih,v)-f(ih,y)) 

D"p h ({i + \)h, kh, v, y)dv for \v\ =3, 
J 5 = h^h / A r -x(0,ih,x,v)(f(ih,v) - f(ih,y)) 
D"ph((i + l)h, kh, v, y)dv for \v\ =4, 

fc-2 „ 

J 6 = h 1/2 J2 h / A r _ 1 (0,ih,x,v)(f(ih,v) - f(ih,y))(kh-ih) 
D» +e ? +e «p(ih,kh,v,y)dv for \v\ = 3. 
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It follows from (O that 



|J 4 | < C^/^Bi 1 -, a) x ... x 5(1 s+ L-l) 



xB(e, e + 



r — 1 



)(fc*0 



C5jgs(y-a:)- 



2 /v ' 'vfeft 
Clearly, the same estimate holds for J5 and J§. Thus we obtain 



)^(e,e+^)(^) 2£+ ^C5^(2/-a:). 



(80) |J/|<C p /i 8 / a - 2e B(ie)x...xB(ie+ 2 
Now we estimate J//. We denote 

Br-i = P ®h (H + M M + VhN^-V -p® h (H + Mh^-V. 
Using the recurrence equation Bq = and 

B r _i = 5 r _ 2 ^ (H + M M + VhNi) +p® h (H + M M ) (r - 2) ® h VhN x 
we obtain that 

r-2 

(81) /// = ®h VhN, ® h (H + M M + VhN.f- 1 ^ ® h M h , 2 (0,T,x,y), 

1=0 

where p 3 j = p®h {H+M^ip 1 '- To estimate III it's enough to estimate a typical term 
in the last sum. Thus we need a bound for 



(82) 



n-2 „ ( fe-1 [" „ j-l „ 

ft3/2 E"/ E" /E"/ 

^=0 ^ L j=o L »=o 



p 3 j(0,ih,x,w)(jh - ih) 



xD^ +en+em p(ih,jh,w, z)(g(ih,w) - g(ih, z))dw\ 

x(H + M M + VhN^-^iih, kh, z, v)dz} (f(kh, v) - f(kh, y)) 

xD»p h {{k + l)h,T,v,y)dv. 

For an estimate of (jEZj) we apply two times integration by parts in the internal integral f ...dw 
and then we make two times an integration by parts in f ...dv. We use also the following 
estimates for < / < r — 3 

(83) \D a w D% tl (0,ih,x,w)\ 

<C^B{l\)x...xB^\)iih) 1 -^ 

(84) D a v D h z {E + M h>1 + y/hN^ r - l - 2 \ih,kh,z,v) 
l B( 



w — X) 



<0<-'-^,i)x...x B (i,^ 



r-l-4-\a\-\b\ 

)(kh - ih) a (Vkh=ih^ 



v — z) 



We put here -B(|, 0) = 1. This gives the following estimate for < / < r — 3, r > 2 
(85) |p 3) i ® h Vhm ® h {H + M M + VhN^ 1 - 1 - 2 ^ ^ M A)2 (0, T, ar, y) 

/2- 3£ r(e) 



r(3 £ + £f± 
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For I = r — 2 we have to estimate 

i?3,r-2 ®h VhN lh <g> h M h<2 (0, T, x, y). 

This is a finite sum of terms corresponding to the different summands in jVj.^ and 
A typical term can be bounded by 



n-2 



' fc-1 



(86) h 3/2 J2 h [ \j2 h [ P3,r-2(0,jh,x,w)(kh-jh)D^ J +e " +em p(jh,kh 1 w,v) 

x {g{jh, w) - g(jh, v))dw} (f(kh, v) - f(kh, y))D v v p h {(k + l)h, T, v, y)dv. 

We apply again integration by parts and after direct calculations we obtain the following 
estimate for r > 2 



(87) 



P3,r-2 ®h ^/hN lh ® h M hj2 (0, T, x, y) 



< C r h 3 / 2 ~ 



3c- 



r(c + ^) 



The inequalities @2||, dH follow from (T70l. (l85l) and (1871). 

Asymptotic treatment of the term T$. We will show that 



(88) 



n 



Y,V®h{H + A) « (0,T,x,y)-J2P ®h # (r) (0, T, x, y) 



r=0 



r=0 



< 



Chn- 5 CMv 



x) 



where A = Mfc — M h = -|(L^ - ILL + L 2 )X(x). Denote 

C r = p^hiH + M'^ + VhN^^T^x^y) 

-p ® h (H + M h + v^iVi) (r) (0, T, x, y) 
-[p ® h (H + A) W - p ® ft ffW] (0, T, x, 2/). 

Analogously to (jB}) we have the recurrence relation 

(89) C r = C r ^ ® h H + \p ® h (H + M'l + VhN^-V 

-p ® fc (H + M h + VhN^-V] ® h [Ml + VhN x ) 



+ 



p ® fe (# + M h + v^AO^- 1 ) — p ® h (H + A) {r -^ 



A 



I + 11 + III. 



Denote 



D r -i = V ®h (H + M h + VhN^-V -p® h (H + Af r ~ l \ 

Clearly 

D r -i = A-2 ®h {H + M h + VhNj + p h ® h (H + Af r - 2) ® h (M h — A + VhN x ) 
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Iterating we obtain 
(90) /// = A_i ®h A 

r-2 



®h (M h -A + VhNi) ® h (H + M h + VhN^- 1 -^ ® h A(0, T, x, y), 



1=0 



where p^i = p®h{H + A)^ . The sum (|9T)|) can be estimated exactly in the same way 
as the sum (JUJ). This gives 



(91) 



\III\ < C{e)ti" 2 ~ 2s ° T " ; 



r(^) 



-T^ + —( VT (v-x),r = 2,3... 



x). 



To estimate II denote 

E r -i = P ®h (H + Ml + VhNxf-V _p® h ( H + Mh + VhN^-V. 
For r = 2 we have E\ = p®h A and analogously to (jlHJ) 
(92) 

Analogously to (IHOIl for r > 3 we use the recurrence relation 

(93) # r _! = E r - 2 ® h (H + Ml + VhN x ) + [p ® h {H + M h + VhN x f- 2) ] ® h A 
= I' + 11'. 

Both terms in l{9~Hj) are analogous to the corresponding terms in (f72T) and may be 
estimated analogously. This gives the following estimates for r > 3 



r-2, 



Er-i\ < C r h 1 - s B(-,e)x...xB(~,s+—-)(khT + ^(^(y 



x) 



(94) 



\II\ 



E r _ x ® h (Ml + VhNj (0, T, x, y) 



< C&h^^T^t^v-x) 



The desired estimate (JHHJ) follows from ((HJ), JM} and (1941) . 

Asymptotic treatment of the term T 4 . We will show that 



(95) 



T 4 = J2P®hH^(0,T,x,y) 

oo 

^ p <g> h [H + /iiV 2 ] « (0, T, x, y ) + it!* (x, y ) , 



r=l 



r=l 



with -/V 2 (s, t, x, y) = (L — L)-k 2 {s, t, x, y), \Rl(x, y)\ < Chn 5 (^{y — x) for 5 > small 
enough and with a constant C depending on 5. For the proof of (|95|) it suffices to show 
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that for 5 small enough 
(96) 



®h (H + M h + Vh^ + hN 2 f'\0, T, x, y) 

r=l 

n 

-J2P®h (K h + M h )^(0,T,x,y) 



< 



E 



r=l 
C k 



hn 6 (%(y-x), 



(97) 



k=l ^2' 
n 

Y J P® h {H + M h + VhNjW (o, T, x, y) 



r=l 



r=l 



^ p ® ft (# + M h + v^iVi + hN 2 ) (r) (0, T, x, y) 
=i 

n n 

P ®hH^ (0, T, x, y) - J2P ®ft [H + (r) (0, T, x, y) 



< 



E 



r=l 



r=l 



k=l ^2> 

Denote D 3 = and 



C/in-*C5r(j/-x). 



D 3>m (Q,jh,x,y) = J2P^h{K h + M h f\0,jh,x,y) 

r=l 

m 

~Y,P®h (H + M h + VhN x + hN 2 f\0,jh,x,y). 

r=l 

Then (|96|) can be rewritten as 

\D 3tn (0,T,x,y)\ < Chn^C^iy - x). 
We now make iterative use of 

(98) D 3>m = D 3 , m _i ® ft (# + M ft + VhNi + hN 2 ) + £ m _ x , 

for m = 1, 2, where 



g m (0,jh,x,y) 



r) 



r=0 



(H-K h + VhN x + hN 2 )(0,jh, 



Sh,m ®h {L - L)d h (0, jh, x, y) 



with 



g (0,jh,x,y) = -p ® h (H - K h + Vh^ + hN 2 )(0, jh, x, y), 
d h = Ph-P- v^tti - hn 2 , 

m 

S h , m (0,ih } x,y) = ^p® h (K h + M h ) {r \(},ih,x,y). 

r=0 
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Iterative application of (|98|) gives 

n-l 

D 3 ,n(0, T, x, y) = 9r ®h (H + M h + VhN, + hN^-^Xo, T, x, y). 



r=0 



To prove (|96|) we will show that 



(99) 



g r (g) h (H + M h + v^iVi + hN 2 ) {n - r - l) (0, T, x, y) 



For this purpose we decompose the left handside of (jlOOj) into four terms 

h J g r (0,ih,x,u)(H + M h + VhN^ hN 2 ) {n - r - 1 \ih,T,u,y)d 

Sh, r (0, kh, x, v)(L — L)dh(kh, ih, v, u) 



a r2 



<V3 




CL r 4 



0<i<n/2 

E * 2 E 

n/2<i<n 0<k<i/2 

x(H + M h + VhNi + hN 2 ) (n - r - 1) {ih ) T, u, y)dvdu, 

/i 2 y / (L T — L T )Sh, r (0, kh, x,v)dh(kh,ih,v,u) 

n/2<i<n i/2<k<i~n s ' 

x(H + M h + \fhN x + hN 2 ) {n - r ~ l \ih, T, u, y)dvdu, 

h 2 y J (L T — L T )Sh, r (0, kh, x,v)d h (kh,ih,v,u) 

n/2<i<n j_ n <S'<fc<j_i 

x (if + Af fc + VhNi + hN 2 ) {n - r - l \ih, T, u, y)dvdu. 



Here L T and L T denote the adjoint operators of L and L, and 5' satisfies inequalities 

K 

show for I = 1,2,3,4 
(100) 

< 



2x < 5' < |(1 — x), where x is defined in (B2). For the proof of ()100p it suffices to 



|ar.i| < hn 5 C n - r B{l, 1 -) x ... x B( U r 2 1 ,^)C%(y-x) 



C n- 



hn-'t^y - x) 



for some 5 > 0. 



Proo/ o/ (EH /or Z = 2. Note that fc < i/2, i > n/2 imply ih - kh > J. The 
claim follows from the inequalities 

(101) max{\K h (ih, jh,x,y)\ , \M h (ih, jh,x,y)\ , Vh~Ni(ih,jh,x,y) 

\hN 2 (ih,jh,x,y)\ , \H(ih,jh,x,y)\} 
< Cp _1 C p (?/ — x) with p 2 = jh — ih for < i < j < n, 
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(102) 
(103) 

(104) 



\S h>m (0,kh,x,v)\ < C( s ^(v - x), 



(L - L)d h {kh, ih, v, u) < Ch 3/2 {ih - khy 2 ( s ^_ kh {u - v) 



0(hri 



-l/2+3/2x\ 



[u — v) 



\Jih~kh 

(H + M h + VhN t + hN 2 Y n - r -V(ih, T, u, y) 

< cr~(f~- 3 B(~, ~) x ... x B( n ~ r 2 ~ 2 , ~)C^(y - u) 



< 



C r, 



n—r—l s 



for n — r — 3 = —1, 0, 1, n — 3 with p 2 = T — ih. We put -B(^, 0) = 1). Inequality 
(jlOlj) follows from the definitions of the functions Kh,...,H. Inequalities (|102j) and 
(110411 can be proved by the same method as used in the proof of Theorem 2.3 in 
Konakov and Mammen (2002) (pp. 282 - 284). Inequality flTTIHl) follows from the 
inequality ih — kh > ?, Lemma 5 and the arguments used in the proof of Lemma 7. 

Proof of JM3) for I = 3. Note that n/2 <i,k> i/2 imply kh > \. We use the 
following inequalities 

(105) \d h (kh,ih,v,u)\ < Ch 3/2 (ih - kh)~ 3/2 ( 



5-6 



(L 1 - L 1 ) V(0, kh, x, v) < CT-\^(v - x) 



\/ih—kh 
1/-5-2 



U — V) 



(106) 



i/2<k<i-n s ' 



h (L T — L T )Sh, r (0,kh,x,v)dh(kh,ih,v,u)di 

(ih-kh) 3 / 2 ' 

x)<Ch 3 l 2 T- 3 l 2 n^ 2 C vm {u-x) 



< C h 3 ^ y: K ih _\ h w CM»-* 



i/2<k<i-n s ' 

" ih - nS ' h du 



^Chn^'t^iu-x), 
where 5" = 5'/2 - x > 0. Claim (fmnjl for I = 3 now follows from (fTO!^ and (fTH^ . 

Proof o/ | f77^) /or / = 4. For i - r/ < /c < i - 1, n/2 < « we have i/i > T/2, 
/c/i > T/3, (i — k) < n s for sufficiently large n. The integral 

(L T — L T )Sh, r (0, kh, x, v)ph(kh, ih, v, u)dv 



is a finite sum of integrals. We show how to estimate a typical term of this sum. The 
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other terms can be estimated analogously. We consider for fixed j, I 

d 2 S htr (0, kh, x, v) 



(107) 



dvjdvi 



(aji(kh,v) - o-ji(kh,u))h d/2 



i-l 



xg (i k) [kh,u,h 1/2 (u - v - h ^ m(lh, u))]di 
d 2 Sh, r {0, kh, x, v) 



l=k 



dvjdvi 



x[(jji(kh,u* — Vhw) — (Tji(kh,u)]q( l k \kh,u,w)dw, 

where u* = u — h Ym=1 m {lh, u). Now using a Tailor expansion we obtain that the ri 
hand side of ()107j) is equal to 



d 2 S htr (0,kh,x,u* 
dvjdvi 



f 1 D u &S h<r (0, kh, x, u* - 5 Vhw) d§ 
v\ J v dvidv) 



UjUVl 



<jji(kh,u) 



M=l 



+2hJ2 



[w + VhJ2izlm(lh,u)Y 



\u\=2 
1 



i-l 







D v u a^(kh, u — 5Vhw — 5h m(lh, u))d5 



l=k 



] {i k \kh,u,w)dw. 



Note that 



-Vh J ° 2 ^ X ' ^ (w p + VhJ2m p (lh, u))q^ k \kh, u, w)dw 

J l=k 



i-l 



d 2 S htr {0,kh,x,u*) ^ 
—h — — > m p {lh, 



(108) h 



dvjdvi 

d 2 Sh, r (0,kh,x,u* 
dvjdvi 



l=k 



i-l 



i-l 



[w p + Vh'^^ m p (lh, u)){w q + Vh"^^m q (lh, u)) 



l=k 



l=k 



D v u <jji{kh, u)d5 + J ^D^aji(kh, u — 5\fhw 



i-l 



—5h m(lh, u)) — D v u Oji{kh, 



u 



l=k 



d5 > q {i k) (kh,u,w)dw 



h d2Sh ^ k ^ ^ Dfaikh, u) J w p w q q^ k \kh,u,w)dw 
+h 2 fe ' r ^ ' n ' — - — - m p (lh, u) m q (lh, u) + R, 



dvjdvi 



l=k 



l=k 
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where by (A3') we have for jo < (i — k) < n 5 ' ,w' — (i — k) l l 2 w 



\R\ < Ch 3/2 



d 2 S h ^ r (0, kh, x, u* 



dvjdvi 



J (>72 +0 {T 1 / 2 n- 1 l 2+s '))%{w')dw' 



(109) < ChC^u - x)(h s n ss ' + i)T- 3 /2 n -i/2+35'/2 J \\ w 'f^ w >) dw > 

< Chn^l^PQ^u -x)< ChnW-^Q^u - x), 
We obtain analogously 

" L d 3 S htr (0, kh, x, u* - dVhw) 



i-i 



h / w p (w q + \fh S ^m q (lh,u))D e ^a j i(kh, 



l=k 



dv p dvjdvi 



d5 



= h 



xg (l k \kh,u,w)dw 
d(jji(kh, u) d 3 Sh, r {0, kh, x, u*) 



dvpdvjdvi 



w p w q q^ k \kh,u,w)dw + R, 



where 



\R\ < Chn-^-^'-^C^u -x),l- 35' - 3x > 
and, for 1 - 35' - 2x > 

%1 d 3 S hjr (0, kh, x, u* - 5y/hw) 



(110) 



h 3 ' 2 



w 



-dS 



dvpdvjdvi 

1 d 2 <jji(kh,u — 5Vhw — ShY^iZl m q(lh-> u )) 



du r du s 



d5 



i-l 



i-1 



x(w r + Vh S ^m r (lh,u))(w s + Vh^2 m s(^h,u))q^ k \kh,u,w)dw 

<Chn-^- 36 '- 2 "\J h {u-x). 
For 1 < i — k < jo the same estimates remain true because the following bound holds 



(111) 



J \\ W \\ S Q^\t, x, w)dw < C(Jq). 



The same estimates hold for p(kh, ih, v, u) with <p^~ k \kh,u,w) instead of q^~ k \kh, u, w), 
where <p(kh, u, w) is a gaussian density with the mean and with the covariance 
matrix equal to a(kh,u). The first two moments of g(*~ fc ) and coinside so after 

substraction we obtain uniformly for % — n 5 ' < k < i — 1 



(112) 




(L T -L T )S h , r (0,kh,x,v) 



E ^ E 

n/2<t<n i—n 5 ' <k<i-l 

x (j> h (kh, ih, v, u) — p(kh, ih, v, u))dv 
x(H + M h + v^iVi + hN 2 ) {n ~ r - 1 \ih,T,u,y)d 



U 



< 



' n—r—1 • 



hT V2 n -m-«-ss>m Cv _( u _ x y 
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To estimate the other terms in dh{kh,ih,v,u) we need bounds for the following 
pressions 

h f{L T -L T )S h>r (0,kh,x,v)Vh{ih- kh) 



i—n s '<k<i—l 



We have 



(113) 



xD%p(kh,ih,v,u)dv for \u\ — 3, 
h Yl f(L T -L T )S h>r (0,kh,x,v)h(ih-kh) 

i—n s '<k<i—l 

xD v v p{kh, ih, v, u)dv for \u\ =4, 
h J {L T -L T )S Kr {Q,kh,x,v)h{ih-khf 

i—n 6 '<k<i—l 

xD l/ u p{kh,ih,v,u)dv for \u\ = 6. 



h J{L T - L T )S h) r{0,kh,x,v)Vh{ih-kh) 

i~n 6 ' <k<i-l 

D v v p(kh, ih, v, u)dv\ 
h Yl j D e ^ +e "{L T -L T )S Kr {Q,kh,x,v)Vh{ih-kh) 



i—n s ' <k<i—l 

D"~ £p ~ eq p(kh, ih, v, u)dv 



h 



-n s <k<i-l 



y/ih — kh 



KChn-^-^C^u-x). 



Clearly, the same estimate (I113j) holds for \u\ = 4 and \u\ = 6. Now (jlOOjl for I 
follows from this remark and (j!12l) and (11131) . 

Proof of nUm for 1 = 1. Note that for this case T - ih > T/2. 
(114) a nl = J2 h " E / f(L T -L T )S h , r (0,kh,x,v) 

0<i<n/2 0<k<i-l J 

xdh{kh, ih, v, u)\^h, r {ih, T, u, y)dvdu 



h J{L T - L T )S h)r (0, kh, x, v) 



0<fc<n/2-l 



x < Y^ h J dh{kh,ih,v,u)^h,r{ih,T,u,y)du 

{k+l<i<k+n s ' 

+ h J d h (kh,ih,v,u)^> h,r{ih,T,u,y)du\ dv, 

k+n s ' <i<n/2 ) 
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where we denote 

* h>r (ih,T,u,y) = (H + M h + VhNi + hN 2 ) ( - n ~ r ~ 1 \ih, T, u, y). 

We consider 



J2 hi h- d ^ 2 q^- k \kh,u,h- 1/2 [u - v - h^2m(lh,u)])V hi r(ih,T,u,y)du 

<i<k+n 5 ' l=k 

J2 hf | q i - l - k \kh,v,w) + VhJ2( w + VhJ2rn(lh,u)yD u v q il ~ k \kh,v, 



IH=2 



l"l=3 



1 i-1 

(1 - 5fD v v q il ~ k \kh, v + <5/i 1/2 u> + <J/i ^ m(Z/i, u), w)d<J 



i=k 

i-1 



X 



1^1=1 z=fc 



M=2 



+3/l 3/ 2 y- ( m + ^Ell m ( f M)f 
^ v\ 

M=3 

x / (l-5) 2 D^ )r (i/i,r,v + 5/i 1/2 w + 5/i^m(Z/i,w),y)^ldw 

Z=fc J 

This integral is a sum of 4 x 4 = 16 integrals. We estimate only two of them. Other 
integrals can be estimated by similar methods. First, we estimate 

h J q^ k) (kh,v,w)^ h:r (ih,T,v,y)dw= ^ h* h , r (ih,T,v,y)dw. 

k+l<i<k+n s ' k+l<i<k+n s ' 

Note that we get the same term when we replace q ( -' l ~ k \kh,v,w) by (p^~ k \kh,v,w). 
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After the replacement this term disappears. Second, we estimate 

i-l 



J2 h j q (l ~ k \kh, v,w)VhJ2( w + VhJ2 m (lh, u)f D U V <S> Kr (ih, T, v, y)dw 

l<i<k+n s ' |f |=1 l=k 

d „ i—1 

^E E D^ r (ih,T,v,y) J q^ k \kh,v,w)[w, + VhJ2 m ^ v ^ 

3 = 1 k4-l<i<k4-n s ' " l=k 



J =1 k+l<i<k+n< 

+0(hn 6 ' \\w\\ + h 3/2 n 26 ')}dw 

d i-l 
i =1 k+l<i<k+n 6 ' l=k 

+o(h 2 n 25 'J2 E h\DZ* h , r (ih,T,v,y 

\ i =1 k+l<i<k+n s ' 

+olh 3 VJ2 E h\Dl'^ htr (ih,T,v,y)\Jq^- k \kh,v,w)\\w\\dw\ 

\ J =1 k+l<i<k+n s ' / 

d i-l 

= h " E E D v j *h,r(ih, T,v,y)J2 mj (lh, v) + R, 

3=1 k+l<i<k+n s ' l=k 

where 

(~in—r 

\R\ < Y^T l/2 hn-^\^ h {y-v). 

The first term in the right hand side of this equation will be the same if we replace 
q^ l ~ k >{kh, v, w) by <p {i ~ k \kh, v,w). After the replacement this term disappears. For 
a proof of this equation we consider the function u(w) that is defined as an implicit 
function and we used the following change of variables 

i-l 

b}l 2 w = u — v — h m(lh, u) 
i=k 

to obtain 

i-l i-l 

\fh^m{lh,u{w)) = VhJ2 m ( lh > v ) + (M* - k ) W W W + h3/2 ^ ~ k f) 

l=k l=k 

because of (i — k) < n 5 ' . By similar methods we get 



(115) 



h [VhTTi(kh,ih,v,u) + h7T2(kh,ih,v,u)]^h,r(ih,T,u,y)du 



k+l<i<k+n 



<^^hn^ 2+2S '^ 2 (^ h (y-v). 
It remains to estimate 

N iu. 



^ hi d h (kh,ih,v,u)^ h . r (ih,T,u,y)di 



k+n s ' <i<n/2 
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From ijTTO and (IT051) we obtain 
(116) 



h / dh{kh,ih,v,u)^h,r{ih,T,u,y)d 

T -l/2 h 3/2 



k+n 5 ' <i<n/2 


< 


(jn—r 


T (n-r-l 


< 


(J n ~ r 


V (n-r-l 


< 


(~jn—r 


T (n-r-l 



T/2 



kh+n s 'h ( u ~ kh) 3 / 2 



Now we substitute the estimate (|116|) into (|114|) . This gives the following estimate for 
any < e < x 



(117) 



J2 h f(L T - L T )S hir (Q,kh,x,v) 
k=i ^ 

h J df l (kh,ih,v,u)tyf ltr (ih,T,u,y)d 



< 



k+n 6 ' <i<n/2 
(jn—r 



re 



n—r—l ' 



n/2 



<C(e) 



r( 



n—r—l 



hn -l/2(5>-x) h -e J- hikhf^C^iy ~ 
k=l 

hn ~l/2(S>-»s) Cv _ {y _ x) . 



X) 



For fc = 0we get with Sh, r (0, 0, x, v) = 5(x — v) where £(•) is the Dirac function that 

Sfe, r (0, 0, x, v)(L - L)d h (0, ih, v, u)^ h ^ r {ih, T, u, y)du 



E " 2 

l<i<j/2 




<C(e) 



' n—r—l ' 



/in -(V2- E ) Cv _( 2/ _ ;r ). 



This completes the proof (jlOOj) for / = 1. The estimate (1371) may be proved by the 
same arguments as were used to prove (JHEJ)- 

Asymptotic treatment of the term T 5 . We will show that, 

oo 

(118) T 5 = -VhJ2*i ®h (H + M M + VhNjWtp, T, x, y) 

r=0 



-h k 2 ® h (0, T, x, y) + R h (x, y), 



r=0 



where \Rh(x,y)\ < Chn 7 C/y {v ~ x ) f° r some 7 > 0. Note that with Sh(s,t,x,y) 
Ylr=i(Kh + M} l p r )(s,t,x,y) the term T 5 can be rewritten as 

T 5 = (p-Ph)(0,T,x,y) + (p-Ph) ®h S h (0,T,x,y). 
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We start by showing that for x < 5 < uniformly for x,y G R 



(119) 



h ^2 (Ph-p)(QJ h ,x,u)S h (jh,T,u,y)d 

l<j<n s 

<0{hn- l W-*-^)C% 2 {y-x) 



u 



for 5 small enough. For the proof of (111911 we will show that uniformly for 1 < j < 
and for x,y G R d 



n 



(120) 



(121) 



p h (0,jh,x,u)S h (jh,T,u,y)du = S h (jh,T,x,y) 
+0[h 1 l 2 T- l l 2 n- 1 ' 2+s + h}' 2 T~ l + n 5 / 2 ^ 2 ]^ 2 (y - x), 

p(0,jh,x,u)S h (jh,T,u,y)du = S h (jh,T,x,y) 
+0[h 1 ' 2 T- 1 / 2 n- 1 ' 2+5 + h x ' 2 T- 1 + n 5 ' 2 ^ 2 }^ 2 (y - x). 



Claim ()119p immediately follows from (|120|) - (|121j) . For the proof we will make use of 
the fact that for all 1 < j < n 5 and for all x,y G R d and \u\ = 1 

(122) \D:S h (jh, T, x, y) \ < C(T - jh) \Y ,.,.(// ~ *)■ 

Claim ()122|) can be shown with the same arguments as the proof of (5.7) in Konakov 
and Mammen (2002). Note that the function $ in that paper has a similar structure 
as Sh- For 1 < j < n s the bound (I122JI immediately implies for a constant C 



(123) 



iD^ijh^x^^KC'T-'C^iy 



x). 



We have ph(0,jh,x,u) = h d ^ 2 q^[0,u,h 1 l 2 {u — x — ^ELo m (^) M ))]- Denote the 



determinant of the Jacobian matrix of u—h^2 J i=Q m(ih,u) by A h . Because of condition 
(A3) and (IT^l) it holds that for 1 < j < 



p h (0, jh, x, u)S h (jh, T, u, y)du 



n 



j h- d/2 q {j) [0, u, h~ 1/2 {u -x-hJ2 ™(ih, u))]S h {jh, T, u, y)d 

r j ~ l 

I q {j) (0,x + h 1/2 w + h S ^m(ih,u{w)),w) |A^| 
J i=0 

Sh(jh, T, x + h x l 2 w + h m(ih, u(w)), y)dw 



i=0 



[ q ^\o,x,w)+o(r d/2 h i / 2 )(\\w\\ + mr 1/2 wm 

+ 0(h 1 ^ 2 T- 1 )C s VT 2 (y - x)(l + h^' 2 \\w\\ s - 2 )(\\w\\ + l)]dw 
S h (jh, T, x, y) + Oih^T-^n- 1 ^ 5 + h^T' 1 + /W 2 ]C^ 2 (j/ - x) 
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with u = u(w) in X^=o m (^ M ) defined by the Inverse Function Theorem with the 
substitution w = h~ x l 2 {u — x — h^2i~Qm(ih,u)). This proves (11201) . The proof of 
(I121j) is the same with obvious modifications. From (|119j) we get that for 5 < 
(where x is defined in (B2)) 

T S = 9 -mO,T,x,y) + H £ /(?-ft)(0,ifc,»,u )& (,T,«, B )*, + i fc (,, y ) 

n s <j<n 

with |fl h (x,3/)| < Oihn-W 1 -"- A6) )C% 2 {y ~ x). We now make use of the expansion of 
Ph — V given in Lemma 5. We have with p = (jh) 1 ^ 2 > h}' 2 n s ^ 2 



(124) 



n r 
h J2 h 3/2 p- 3 / Cp (U ~ x)S h (jh, T, it, y)efo 

< Ch 2 T-*'n- 5 " P~ 2+25 ' / |C?(«-«W».r,«,y)|dtt, 



where 5' < - 5)' 1 , 2 5" = 5 + 255' - 25'. Now we get that 

n „ 

(125) hJ2p~ 2+2S ' / \C S p (u-x)S h (jh,T,u,y)\du<CB(5',l/2)T 5 '- 1 / 2 C% 2 (y-x) 

j=n" 

for a constant C. This shows that for 5' > small enough 

T 5 = -[^i + ^ 2 ](0,T,x,y) 

-h ^2 / [V^tti + /i7T2](0, jh,x,u)S h (jh,T,u,y)du + R' h (x,y), 



n°<j<n 



with l-R^x, 2/) | < 0(/m *^)C^p 2 (y — x) with a constant in O(-) depending on 5'. 



It follows from ifTTTjj) . (112411 and (11251) that 

oo 

(126) T 5 = - ^[VfcSri + /wr 2 ] ®fc (#h + M h f\0,T,x,y) + R'^x,y), 



r=0 



where {R'^x, y)\ < 0(hn ^" ^ 2 ^)q'^ 2 (y — x). Now we apply Lemma 10 with A 
Vhn u B = H + M hA + VhN u C={K h -H- VhN,) + (M h - M M ) to 



(127) - ^ V^tti ® h (K h + M h ) W(0, T, x, y) 

r=0 

oo 

+ ^ V^tt! ®„ (if + M M + V^iVi) (r) (0, T, x, y) 

r=0 

and with A = hn 2 , B = H,C = (K h - H) + M h to 

oo oo 

(128) - h ^ ®h (K h + M fc )W(0, T, x, y) + ^ /wr 2 ® ft # (r) (0, T, x, y). 



r=0 



r=0 
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The estimate (fTT8j) follows from (fT25l . (fHTI) . (IT281) . Lemma 10, Lemma 5 . ([THl) and 
(HH). 

Asymptotic treatment of the term T§. By Lemma 9 

iTelKC&hn-^C^iy-x). 

Asymptotic treatment of the term T-j. We use the recurrence relation for r = 2, 3, ... 

Ph ®h (K h + M h + R h )W (o, T, ?/) - p ft ® ft (0, T, x, y) 

= [p h ® h (K h + M h + R h y r ~V - p h ® h H^-V] ® h H h (0, T, x, y) 
+ [p h ® h {K h + M h + R h f-V ® h {K h + M h + R h - H h )){0, T, x, y) 
and we apply Lemma 8 for r = 1. We get in the same way as in the proof of Lemma 9 

\T 7 \ < Ch^T" 1 / 2 ^ -x) = ChrT 1 ' 2 ^ - x). 

Plugging in the asymptotic expansions of T 1? ...,T 7 . We now plug the asymptotic 
expansions of T 1; ...,T 7 into (|38| . Using Lemma 10, Theorem 2.1 in Konakov and 
Mammen (2002) we get 

(129) p h (0,T,x,y)-p(0,T,x,y) 

= Vh [tti + p d ® h 3«i] ® h $(0, T, x, y) 

+h { [tt 2 + tti ® h $ <g) h sRi + p d ® fe 3? 2 + / ®fc 32s] ®h $(0, T, x, y) 
+p d ® h (Ki ®,$) (2) (0,T,x,y) 

+ ip ® h (Lj _ L 2 )/ (0, T, x, y) - ~p ® h (L' - L')p d (0, T, x, y) 
+0(h 1+s ( VT (y-x)) } 

where 

9?i(s,f,x,y) = JVi(s,£,a:,y) + M l {s,t,x,y) - M x {s,t,x,y), 
9£ 2 (s, *, s > y) = N 2 (s, a?, y) + ni(s, t, x, y) - fii(s, t, x, y), 
$t 3 {s,t,x,y) = - D xP (s,t,x,y), 

\u\=4 V - 

Mi(s,t,x,y) = Xu(yS, } ^ D xP( s > t, x, y), Mi (a, t, x, y) = ^ ^y^ D xP( s , t, x, y), 

\u\=3 U ' \v\=3 ' 

Ux(s,t,x,y) = V Xu(yS ' ^ D^rri(s, t, x, y), Ui(s, t, x, y) = V Xv S ' V D v x it x (s,t,x,y). 
z — ' v\ z — ' v\ 

\u\=3 M=3 

For the homogenous case and T = [0, 1] (|129|) coinsides with formula (53) on page 
623 in Konakov and Mammen (2005). 
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Asymptotic replacement of p d by p. It follows from (11), (26) and (27) in Konakov 
(2006) that 



\(p d -p)(ih,jh,x,z)\ < C{e)h 1 - £ {jh-ih) £ - 1/2 (j) 



[z — x) 



(130) \yp —p)\ub,jib,u,,A)\^^ye.)ib yjib — ub) v ^J{j-i)h y 

for any < e < 1/2. Using ()130ft and making an integration by parts we can replace 
p d by p in (|129jl . For example the operator L 2 — L 2 is the operator of the third order. 
Making an integration by parts we have for \u\ =3 



(131) 



n— 1 „ 

Y^h I D» z p(0,th,x,z)(p d -p)(ih,T,z,y)dz 
i=i ^ 



i=l 



ihf 2 {T-ihy/ 2 - 



:<t>vr(y- x ) 



< C&h^T^Bie, e + -)^(y - x). 



By (B2) we have < x < 1 — As. This implies 



-p® h (L 2 -L 2 )(p d -p)(0,T,x,y) 



< C(e)hTy A n 



l/2_-(l/2-2e-*/2) 



for some < 5 < 1/2. The other terms in (I129B containing p d can be estimated 
analogously. Thus we get the following representation 

(132) p h (0, T, x, y) - p(0, T, x, y) 

= v 7 ^ + p ® A &i] ® ft $(0, T, x, y) 

+/i {[ff 2 + 7Ti 0ft $ <g) h 3?i + p ® h ^2 + V ®h ®fc $(0, T, x, y) 
+P ® ft ® ft $) (2) (0,T,x,y) 

+ \p ®h (Ll - L 2 )p(0, T, x, y) - ip ® ft (L' - L')p(0, T, x, y) 
+0(/ i ^Qr(2/-*))- 

For the further analysis we need a generalization of a binary operation Cg> introduced 
in Konakov and Mammen (2005). Recall a corresponding definition. Suppose that 
s G [0,t — h] and i G {h,2h, ...,T}. Then the binary type operation ®' h is defined as 
follows 

fK^v)= E *//(.,iM,<) 9 UM, 

s<jh<t~h 

Note that for s G {0, fa, 2/i, T} these two operations coincide, that is ®' h = ®h- 

Asymptotic replacement of {p®h^-i) ®k$(0, T, x, y) £>y p® (3?i<8>/ l $)(0, T, x, y) = 
(p®3?i)®/ l $(0,T,x,y),z = 1,2,3, [p® h (&i ® A $)] ® ft (&i ® h $) (0,T,x,y) 6y p <g> 
Pi^$)®U^S4$)](0,T,x,y),p®,(L 2 -L 2 )p(0,T,x,y) 5y p® (L 2 -L 2 )p(0, T, x, y) 
and p® h {L' -L')p(0,T,x,y) by p®(L' — L')p(0, T, x, y). 
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These replacements follow from the definitions of 3^, i = 1,2,3, and can be proved 
by the same method as in Konakov (2006), pp. 9-12, where we estimate the replace- 
ment error of p®hH by p®H. Linearity of the operation <S>h implies that it is enough 
to consider the functions p <&h S where £s(u, t, z, v) is a function of one of the followin 
forms: 

Xu(u,z) ~Xu(u,v) . Xv{u,z) ~Xu(u,v) Alio 

- y D x p(u, t, z, v), \v\ = 3, 4, — D^u, t,z,v),\u\=3 

(L - L)tti(u, t, z, v),(L — L)n 2 (u, t, z, v). 

We consider the case $f(u, t,z, v) = (L — L)iri(u, t, z, v). The other cases can be treated 
similarly. It is enough to consider a typical term of (L — L)ni{u, t, z, v ), namely, we 
shall estimate 

/jh r / rjh \ _ 

du p(0,u,x,z)i / Xu(w,v)dw J D V Z (L — L)p(u, jh, z,v)dz 

~^^h J p(0,ih,x, z) (^j^ Xv(u;,v)dw^j D U Z (L — L)p{ih, jh, z,v)dz 

i- 1 p(i+i)h p J- 1 p(i+i)h p 

du [A(u) - X(ih)]dz = J2 . ^ u ~ ih "> du / x '( ih ) dz 



1 J -1 p(i+l)h pi p 

(133) +oJ2 {u-ihfdu / (1-5) X"(s)\ s=Si dzd5du, 

2 - =0 Jih Jo J 

where X(u) = p(0, u, x, z) ^ f^ h Xu(w, v)div s j D u z H(u, jh, z, v), = ih + 5(u — ih). As in 
Konakov (2006), p. 9, we obtain that 

J -1 p(i+l)h p 

/ (u — ih)du j X'(ih)dz 

. _n J ih J 



i=0 
rjh 



h / Xu( s ,v)ds / p(0,ih,x, z)D"A Q (ih, jh, z,v)dz 

h J_1 P h f 
+ / Xu{s,v)ds j p(0,ih,x,z)D u z H 1 (ih,jh,z,v)dz 



(134) -- / p(0,ih,x,z)D u z H(ih,jh,z,v)dz, 

1 i=o J 



where 



A (s,jh,z,v) = (L 2 - 2LL + L 2 )p(s,jh, z, v), 
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Hi(s,t,z,v) 



d l a i:j (s,z) d l aij(s,v)\ d 2 p(s,t, z,v) 



d / Ql „ / ,. - i ;)', 



ds l 



ds l 



dzidzj 



d l mi(s,z) d l mi(s,v)\ dp(s,t, z,v) 



ds l 



ds l 



dzi 



-,l = 0,l,2,H = H. 



The differential operator A was calculated in Konakov (2006), p. 9. A is a forth 
order differential operator. From the structure of this operator and from (I134j) it is 
clear that it 'is enough to estimate 
(135) 

h n ~ 1 f 1 P h f 
1 = 7^^2 h / zZ h / Xu(s,v)ds / p(0,ih,x,z)D" +IM p(ih,jh,z,v)dz$(jh,T,v,y)dv 

Z j=0 J i=0 Jih J 

for \v\ = 3, = 3. To estimate (|135j) we consider three possible cases: a) jh > 
T/2, ih < jh/2 jh -ih> T/4 b) jh > T/2, ih > jh/2 ih > T/4 c) jh < T/2 
T — jh > T/2. In the case a) we make an integration by parts transferring two 
derivatives to p(0, ih, x, z). This gives 

[\vMds [ D^ P (0,ih,x,z)D^- e ^p(ih,jh,z,v)dz 

1 i=0 Jih J 



< Ch 1 -^ 



jh 



du 



u 1 -* (jh-u) 1 -^^ 
<C{e)h 1 -*(jh)*- 1 <l> V3S (v-x) 



v — X) 



and 



|/| < C{e)h 



l-2e 



du 

U l ~ 2£ {T - U 



(136) < C(£)/i 3 / 4 T 1 / 2 n-( 1 / 4 - 2£ —/ 2 )0 vT (y -x)< C{e)T 1 l 2 - 5 h i l A+s <P VT {y - x), 

where 5 = (1/4 - 2s - x/2) > if k < 1/2 - As, < e < 0, 05 (see the condition 
(B2)). In the case b) we make an integration by parts transferring four derivatives to 
p(0, ih, x, z). This gives the same estimate (j!36)l . At last, in the case c) we make an 
integration by parts transferring three derivatives to §(jh,T,v,y) and one derivative 
to p(0, ih, x, z). This gives the same estimate ()136j) . To pass from D^p{ih, jh, z,v) to 
D^p{ih, jh, z,v) we used the following estimate 

\D%p(ih, jh, z, v) + D%p(ih, jh, z, v) | < Cfy^^y - z). 

Clearly, the same estimate (I136j) hold true for the other summands in the right hand 
side of 



2 \ r £ i pjh p 

h / / ] h / Xu(s, v)ds j p(0,ih, x, z)D"Hi(ih, jh, z,v)dz$(jh,T,v,y)dv 
j=0 J i=0 J ih J 
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h j Xu(ih,v)22h / p(0,ih,x,z)D u z H(ih,jh,z,v)dz$(jh,T,v,y)dv 



i=0 " i=0 



< C{s)T x l 2 - &1 

Now we shall estimate the second summand in the right hand side of (|134|) . Analogously 
to (19) in Konakov (2006) we obtain 

i i~ l r(i+l)h rl r 

(137) -J2 / {u-ihfdu / (1-5) \"{s) \ s=St dzdSdu 

^ i=Q Jih JO J 

i J -1 r(i+i)h rl 4 rjh 

= oE / ( M -^) 2 / / X.(r,t;)dT 

z i=o ^ - 70 fe=l Js 

x J p(0, s, x, z)D^Ak(s, jh, z,v) \ s=Si dzdddu 

i~ l p(i+l)h rl 



x y p(0, s, x, z).D^4o(s, j/i, z, f ) | s=Sj dzdddu 

3~ l r(i+l)h rl 
,•_„ Jih JO 



x y p(0,s,x, z)D"Hx(s,jh, z,v) \ s=s . dzdddu 

-M> /"n n 8 X-( s .") 
"Si <«-*>./„ A" fl-fiT" 

x J p(0,s,x,z)D"H(s,jh,z,v) \ s=s . dzdddu, 

where the operators Ai,i = 1,2,3,4, are defined as follows: 

A x {s,jh,z,v) = {L 3 — 3L 2 L + 3LL 2 — L 3 )p(s,jh, z, v), 

A 2 {s,jh,z,v) = {LxH + 2LH 1 )(s,jh, z, v), 
A 3 (s,jh,z,v) = [(L - L)L 1 + 2(Li - Lx)L\p(s,jh,z,v), 

(138) A 4 (s,jh,z,v) = H 2 (s,jh,z,v). 

The operator A± is given in (24) in Konakov (2006). As in Konakov (2006) it is enough 
to estimate for fixed p, q, r, I 

j- 1 r (i+l)h r i _4_ r jh 

Xv(r,v)dr 



1 J - r(i+i)n ri i p 

±W in-ihf (1-*)W 
^ i=o ^ ^ fc=l ^ 



pfO s x I ,,., M( /„ 
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As in (25) in Konakov (2006) we obtain that (|139|) does not exceed 



(139) C^h^W^^v-x). 

It follows from the explicit form of these operators (|138j) that the same estimate (|139j) 
holds for A 2 , A 3 and A4. The other three terms in the right hand side of (|137|) do not 
contain the factor p h Xv( T ,v)d,T and they should be estimated separately. Clearly, 
it's enough to estimate the term containing Aq. The remaining two summands are less 
singular. From the explicit form of A (formula (15) in Konakov (2006)) we obtain 
that it is enough to estimate for fixed q, I, r 

(140) Yl / {n-ihf / (l-S) X u(s,v) 

J p(0,s,x,z)D u z ^ - ^Q^dzidl^^ ( s, i h,z,v } ' s=s ' dzdSdu - 
Analogously to (25) in Konakov (2006) we get that (I14()jl does not exceed 

(141) C&h^ijh^-^iv-x). 
Now from ifHEjjl . (fT36l . ljT37j» . (fT39l and (UTTl) we obtain that 

[p ® h (L - Z)tti] <g) ft $(0, T,x,y)-p® [(L - Z)tti <8>' ft $] (0, T, x, y) 



(142) <Ch 3 ^(f>^(y-x) 

for some 5 > 0. The other replacements can be shown analogously. Thus we come to 
the following representation 

p h (0,T,x,y) -p(0,T,x,y) 

= Vh [tt! ®' h $(0, T, x, y) + p <g> (3?! ®' h $) (0, T, x, y)] 
+/i [5r 2 ®' h $(0, T, x, y) + p ® ®' h $) (0, T, x, y) + p ® fe (& 3 $ ) (0, T, x, 1/)] 
+/i [tFj ®' h $ + p ® (3?! ®'„ $)] ®; (^ ®' h $)(0, T, x, y) 

(143) +^p ® (L* - L 2 )p(0, T, x, y)-^p® (L' - L')p(0, T, x, y) + 0(h 1+s {^{y - x)). 
Now we furhter simplify our expansion of ph —p. We now show the following expansion 

p h (0,T,x,y) -p(0,T,x,y) 

= Vh(p ® ^i[p A ])(0, T, x, y) + /i (p ® J- 2 [p A ]) (0, T, x, y) 
+/i (p ® T x \p ® ^1 [p A ]]) (0, T, x, y) 
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(144) +^p ® (Ll - L 2 )p(0, T, x, y)-^p® (L' - L')p(0, T, x, y) + 0(h l+s CMv - *))> 
where for s e [0, t - h], t e {h, 2h, T} 

p A (s,t,z,y) = (p®' h $)(s,t,z,y) 

= p(s,t,z,y)+ h J P( s J h ^ z ^ v )^iU h ^^ v ^v) dv - 

s<jh<t—h 

Here $i = if + H ®' h H + H ®' h H ®' h H + ....We start from the calculation of p g 
L7Ti(s,i,x,y). 

p®Lit 1 (s,t,x,y)= I dr ( p(s,r,x,v)(t - r) ^ Xi/ T '/' y D v v (L v p(T,t,v,y))dv 



l"l=3 



-^2~i[ dv [ P( s , T i x i v )(f Xu(u,y)du\ ^-D^p(T,t,v,y)dn 



(145) 



3 + t ^ 

|z/|=3 7 >/s M=3 7 



Integrating by parts w.r.t. time variable we obtain 



[ dv P( s i r i x i v )([ Xu{u,y)du\ D^p(r,t,v,y) \ T r : 
M=3 l/ ' 7 L V>/t ' 



=(s+t)/2 



a + i 
2 



Op^Sj T, X, 'u) , 

^X^, u, y) ( — / y)^ - 7", x, u)x„(r, y) ) dr 



\u\=3 J 



p( s > ^y^> u) Xi/(«> y)du j L>>(^-Jp, t, y) 



M=3 



|i/|=3 Z/ ' Js VJr ' 



(146) 

x 



/ ^, r, x, .)Z»>, t , „, ,)*,- E I r X „( T , „)* r, x , v)D"J>(t, t, v, ,)*,. 

|^|=3 Z/ ' ^ s ^ 

Analogously we get 

77 =X!~7 / p ( s > ( / t Xu(u,y)du j D"p(?-^-t,t,v,y)dv 

\v\=3 V ' J ~ J 
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+ ^Z^][ s+t dT \[ Xu(u,y)duj I L T p(s,T,x,v)D"p(T,t,v,y)dv 

|i/|=3 V ' ' ^ 

(147) - ~[ / Xu(r,y)dr J p(s,T,x,v)D„p(T,t,v,y)dv. 
From (IH5T) - (fTITIl we have 

(148) LiTi(s,t,x,y) = 7Ti(s,t,x,y) +p® L?ci(s,t,x,y) —p® Mi(s,t,x,y). 
and from (114811 we obtain 

tti(s, t,x,y)+p® 3£i(s, t, a;, y) 
= 7ri (s, t, x, y) + p ® Lttx (s, t,x,y) — p® Ltti (s, t,x,y) + p® M 1 (s, t, x, y) 



(149) 



-p ® Mx(s, t, x, y) =p® Mi(s, t, x, y). 



It follows from (|149jl and the definitions of the operations ® and Cg/^that (x[ s ; j^] 
below denotes the indicator of the interval [s,jh]) 



(150) 



Vh [tti <g>' h $(s, t, x, y) + (p <g> ®' A $(s, *, x, y)\ 

= Vhfa + p ® 3?i) ®' fc $(s, t, x, y) = >//i(p ® Mi) ® ft $(s, t, x, y) 

= Vh h j{p® M 1 )(s,jh,x,z)<&(jh,t,z,y)dz 

0<jh<t-h J 

= Vh hi I du I p(s,u,x,v) 

0<jh<t~h J *-Js J 

xM%(u, jh, v, z)dv] $(jh, t, z, y)dz 

= Vh j2 h 

0<jh<t-h 

xMi(u, jh, v, z)dv] $(jh, t, z, y)dz 




dux[s,jh] / p(s,u,x,v) 



Vh du p(s,u,x,v) 

^ s J \v\=Z 



dv 



h x[s,jh] / p{%jh,v,z)$(jh,t,z,y)dz 

,0<jh<t-h 

sTk[ du J P (s,u, X ,v) 

x ^ D vPA(u 7 t 7 v,y)dv = Vh(p® J"i) [p A ](s, t, x, y). 



M=3 



Using similar arguments as in the proof of (11501) one can show that 
h [tt 2 ®' h ®(s, t, x, y) + (p® to 2 ) ®' h ®{s, t, x, y) 
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(151) +p ® h (3? 3 ®' h $) (s, t, x, y)\ = h{p ® F 2 ) [Pa] (s, t, x, y) + hp® III ®' h $(s, t, x, y). 
For the first two terms in the right hand side of (|143|) we obtain from (|150|) and (|151|) 

Vh [nt ®' h $(0, T, x, y) + (p <g> 3?i) ®' fc $(0, T, x, y)\ 
+h [n 2 ®' h $(0, T,x,y)+p® (& 2 ®' h $) (0, T, x, y) + p ® fe (& 3 «4 $ ) (0, T, x, y)] 

(152) = V% ® ^i) [p A ] (0, T, x, y) + fc(p ® JS) [pa] (s, f , ar, y)+hp® IIj ®' h $(s, t, x, y). 
Using (I15()j) we also have 

/i [tti <g>' ft $ + p ® (Ki ®; $ )] ®' h ®; $xo, t, x , y ) 

= h(p® ^[p A ]) ®' h (3?! $)(0, T, x, y) 
= hp®F 1 [p A ®; (5fti ®' h $)] (0, T, x, y). 

Note that 

fcp®ni®^(s,t, a;, |/) = /» / rfw f p(s,u,x,v)^^^D^ 1 ®' h <S>](u,t,v,y) 

Js J \u\=3 V - 

= hp®?^! ®' h <£>](s,t,x,y). 
For the proof of (|144|) it remains to show that 

hp ® T x [ttl $ + p A (&i ®' h $)] (0, T, x, y) 



(153) = (p ® ^[p ® ^i[p A ]]) (0, T, x, y) + 0(/i 1+ ^(y - x)). 
We will show that 

(154) hp ® ^[(p - p A ) ®' h (3?! ®' fc $)](0, T, x, y) = 0(h 1+s CMv " 
and 

frp®^®'^! ®;$)] (0,T,x,y) 



(155) -/ip ® J^[p ® (s^ ®' h $)] (0, T, x, y) = O^Qr^ - a;))- 

Then (fTTvfl) will follow from (fT54]l . ifTKKl) and flTBTl . The estimate (fTKKll can be shown 
analogously to ()142p . An additional singularity arising from the derivatives in the 
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operator is neglected by the factor h in (|155j) . To estimate 1)1 54J) note that from 
the definition of 3?i and $ 

(156) KaW^XAT,^)! <C{e)h- e {T-jh) e - 1 <i> y/T = 3K . 

Then we use the following estimate which can be proved by the same method as in 
Konakov (2006), pp. 8-12, where an estimate for (p — p d )(0, jh, x, y) was obtained. 

(157) \{p - PA )(u,jh, v,z)\< Ch^tj^iz - v). 
From (ITKHl) and (flSTI) 

(158) \{p- p A ) ®' h (^ ®' h $) (it, T,v,2/)| < C(£)/t 1/2 ~ £ (T - n) £ 0^(y - «) 



Now to estimate 1)154)1 it's enough to estimate a typical term of a corresponding sum, 
namely, for \v\ = 3 we have to estimate 

hi du I p{0,u,x,v) Xu ^ V ^ D u v { ^ h [ (P~PA)(u,jh,v,z) 

J ° J "' {j:u<jh<T-h} J 

,T/2 pT 

x{^ 1 ®' h ^){jh,T,z,y)dz}dv = h ... + h ... = 1 + 11. 

JO </T/2 

To estimate 7J we make an integration by parts transferring three derivatives to 
p(0, u, x,v) Xu ^\' v . Using 1)158)) we obtain the following estimate 



II\ < C{e)0 2 -* jT ^j^du^y - x) 



(159) < C(e)h 3 / 2 - £ T s (j)sr(y - x). 

To estimate 7 we consider two cases, namely, a) jh — u > T/4 and b) jh — u < T/4 
=>- T — jh > T/4 . Analogously to (11) in Konakov (2006) the difference /t(p — Pa) 
can be represented as 

h(p - p A )(u,jh,v, z) = h{p® H — p®' h H)(u,jh,v,z) 

+h(p ® H - p ®' h H) ®' h ^ x {u,jh,v,z) 

— h I dr j p(u,T,v,z')H(r,jh,z',z)dz' 

i- 1 r(i+l)h 

£ 



j- 1 p(i+X)h t- 

+h dr (K T , z ') ~ K ih , z ')) dz + h{p®H-p ®' h H) 

„•_„•* J ih J 



(160) = /' + W + III'. 
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where A(r, z') = p(u,r,v, z')H(t, jh, z' , z),$i(ih, jh, z'z) = H(ih, jh, z'z) + H®' h 
H(ih,jh,z'z) + = j*(u) = [|] + 1 (with a convention [x] = x — 1 for integer 

x). For J', case a), we have jh — r > T/5 for n large enough. With a substitution 
v + v' = z' we obtain 

rj*h r 

D"h / dr p(u,T,v, z')H(t, jh, z' , z)dz' 
ri*h p 

D v v h \ dr p(u,r,v,v + v')H(t, jh,v + v' , z)dv' 
j * h dr 



< Ch 



/jk=u(z-v) <ChT <f>^=z{z-v) 



:i6D 



< Cn (fr^jh^iz -v) = CT h (p^jj^^z - v) 



To get ()161j) we used an estimate from Friedman (1964) (Theorem 7, page 260) 



\D v v p{u, r,v,v + v') I < C(t - u)- d/2 exp 



C\v'\ 



T — U 



For J I'(u,jh,v, z)(3t 1 ®' h &)(jh,T, z,y)dz, case b), it is enough to estimate for \v\ — 3 



t, [ J * h a [\ ( . i\( ( , >\ ( uiff d 2 p(r,jh,v + v',z) , 
hi dr [p(u,r,v,v + v )(cr Zfe (r,w + v ) - a [k {r, z))\D v/ - , dv 

J u J 



dv[dv' k 



(162) 



x(3?! ®' h $){jh,T,z,y)dz. 



Transferring five derivatives from p to (3?i ®^ $>)(jh,T, z,y) and using the following 
estimate for = 5 

l- D ^'P r ( T > w + *A *) + D zP( T , jK v + v', z ) I 
< C (jh - r)- d l^^ T 
we obtain that (I162j) does not exceed 



z — v — v 



C(j*h - r^T-^^y -v)< ChWT^'n-V-'-W^^y - v) 



(163) =o(h 1+s T 1 - s )<j )VT ^(y-v). 

We used that for any < 5 < 1 x < (see condition (B2)). To estimate 

/ II'(u,jh,v,z)(JR,i ®' h $)(jh,T, z,y)dz we use a decomposition (12) in Konakov 
(2006). The estimate for terms in IV containing the first derivatives X'(ih,z') we use 
identity (14) from Konakov (2006) and similar arguments to already used in estimation 
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of f I'(u,jh, v, ®' h $)(jh,T, z,y)dz. The estimate for terms in II' containing 

second derivatives \"(ih,z') follows from (22) and (23) in Konakov (2006). At last 
for II I' the same estimates hold because of smoothing properties of the convolution 
...® h $i(u, jh, v , z). This implies (|154jl and, hence, the expansion ()144j) . 

Asymptotic replacement of pa by p. We shall compare hp <8> ^[paKO, T, x, y) and 
hp®T 2 \p]{Q,T,x,y) . Note that for 2x < 5 < §, \u\ = 4 



/i / c£u / p(0, u, x, z)xu{u, z)D^p(u, T, z, y)dz 



<ch 1+s (T-h d )-'^My- x ) 



5\-2. 



(164) 



< Ch 



1+5 



11 



2>c 



(Tn* - n^h 5 )' 



and, analogously, 
(165) 



hi du I D"\p(0,u,x, z)xu{u, z)]p(u,T, z,y)dz 

T-h s 



< Ch^-^T^^riy - 



The same estimates hold for Pa(u, T, z, y). Hence, it is enough to consider u G [h s , T 
h 6 }. We consider 

/T~h s p 
du J p(0,u, x, z)xv(u, z)D u z (p - p A )(u,T, z,y)dz 



(166) 

By CEEZ) 



pT/2 i-T-h 5 

h ... + h ... = 1 + 11. 

Jh s JT/2 



\II\ 



T-h 5 



T/2 



du J D z\p(°, u , x , z )xA u i z ))(p - p A ){u,T,z,y)dz 



< Ch^n^iy -x) = Ch^-T^iy 



x 



(167) 

For u G [h 5 , T/2] 



Ch 1+ ^ VT (y-x), 1 >0 



T/2 

h I du I D"\p(0,u,x, z)xv(u, z)](p - pa){u,T, z,y)dz 



(168) 



< C0 2 ~ 5 ^{y - 



x) 
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where condition S < § implies 3/2 — S > 1. It follows from (|164p ~ (|168j) that 
(169) hp <g> ^ 2 [pa](0, T, x, y)-hp® ^ 2 [p](0, T, x, y) = OQi 1 ^ ^{y - x)). 
To prove that 

Tip ® ^i[p ® ^i[p - pa]] = Oih^^iy - x)) 

we consider a typical summand for fixed = 3, 
(170) 

rT 



/i / du p(0, u, x, z)xv(u, z)D v z 



dr / p(u,t,z,v)xv(t,v)D%(p - p A )(r,T,v,y)dv 



As before it's enough to consider u 6 [h s ,T — h s }. The integral in (117011 is a sum of 
four integrals 

fT/2 p f{T+u)/2 

h = h du ...D u / dr 



u 



fT/2 r rT 

I 2 = h du ...D v z / dr 

Jh s J J(T+u)/2 

fT-h s r f{T+u)/2 

h = h I du I ...D v z / dr 

JT/2 J Ju 

pT~h s p rT 

(171) h = h du ...D u z / dr 

JT/2 J J(T+u)/2 

Note that r — u > T/4 in the integrand in J 2 . Integrating by parts w.r.t. v we 
obtain by (fT57) 

(172) \I 2 \<C0 2 -^T^{y-x)). 

Furthermore, u > T/2, r — u > h s /2, T — u > h 5 in the integrand in I4. Integrating by 
parts w.r.t. z we obtain 

(173) \h\ < Ch^ 5 T- l ' 2 (P^(y - x)) < CT^ 2 h 3 / 2 ~^ 2 ~ 5 ^{y - x)), 

where, by our choice of 5, 3/2 — x/2 — 5 > 1. To estimate I 3 we use the representation 
(p-PA)(r,T,v,y) = {p ® H — p ®' h H) (r, T, v, y) 
+ (p®H-p ®' h H) ®' h $i(t, T, v, y) 

ds I p(r, s ) v 1 w)H(s 1 T ) w 1 y)dw 
+ t ^[p ( g ) ' h (H 1 + A )](T,T,v,y) 

Y n ~ l f(i+l)h r-1 4 f 

+ (* - iK ? / C 1 ~ 7) Y) / s ' u ' T ' w > ?/) 

^ i=j» Jih Jo k=l J 
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ih+i(t-ih) dwd^dt, 



(174) 



+(p®H-p ®' h H) ®' h $i(r, T, v, y) 



where j* = j*(r) = [r/h] + 1 (with a convention [x] — x — 1 for integer x) , Hi and 
Ak, k = 0, 1, 2, 3, 4, are defined in Konakov (2006) and 

$i(i/i, i'h, z, z) = H{ih, i'h, z, z) + H ®' h H(ih, i'h, z, z) + ... 

To estimate D v v {j> — p^){j, T, v, y) we note that 

h 



D v v J ds J p(r, s,v,w)H(s,T,w,y)dw 

ds I p(r,s,v,v + w')H(s,T,v + w',y)dw' 



rf*h j 

< ( 'I' I ^ -<t>JT^{y ~v)< Ch 2 - 2S <t>^{y - v). 



r (T-S) 



(175) 
Furthermore, 



h 2 f 

~^ D v h P( T dh,v,w)Hi(jh,T,w,y)dw 



r<jh<T-h 



< 



h 2 



h J D v v [p(T,jh,v } v + w')Hi(jh,T,v + w 1 \y)\dw' 



-<jh<T-h s /2 



+ - 



d „ 

h D u J ei+ek [p(r, jh, v,v + w'))p(jh, T,v + w', y)dw' 

i,k=lT-h s /2<jh<T-h 



(176) 



< Ch 2-2 S(Pv _ {y _ v)+ Ch 2- 5 S/ 2(f) ^_ {y _ v y 



Because of a structure of the operator A (see Konakov (2006)) it is enough to estimate 
for fixed i, I, k 

h7 M h 2 ^ f d 2 p{jh,T,v + w',y) , 

(177) —D v ^ h D w \p(r,jh,v,v + w) dw 

T<jh<T-h J 1 1 

With the same decomposition as in ()176j) we obtain that (|177|) does not exceed 

(178) Ch'-^tP^iy-v). 
From flTTfil) and (ITTSl) we obtain that 



(179) 



y \Dl[p ®' h {H X + Ao)] (r, T, v, y) | < Ch 1+ ^^(y - v) 
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for some 7 > 0. It remains to estimate the last summand in (|174jl . It follows from the 
structure of the operators Ak, k = 1,2,3,4, (see Konakov (2006)) that it is enough to 
estimate 



(!80) ^E/ (* ~ ihf / (1-7) 

2 i=j , Jih JO 



k=l 



dw^dwldw'pdw' j 



Hh+j(t-ih) dw'd^dt 



for fixed i,j,p,q. As in Konakov (2006), p. 12, we obtain that (|180jl does not exceed 

,1 „i n-l 



r r . . 

Ch^^y-v) / z 2 (1- 7 )X> 

JO JO „•_„•* 



1 



[(ih — r) + '-fhz] 3 / 2 [(n — / yz)h — ih} 2 
Ch^^y-v) fz 2 A1-7) £ 



d'ydz 



1=] 

1 



{i:j*/i<i/i<r+/i' 5 /4} 



(181) 



+Ch 2 ^(y-v) I z 2 I (I-7) 



E 



... = I" + 11". 



{i:r+h s /4<ih<T-h} 



1 



< Ch^h-^^iy-v) f z 2 ^ /* ^— ^ 

Jo Jo 7 



{i:j*/i<ih<r+h' s /4} 

E " 

{j:i*/i<i/i<r+h 4 /4} 



[(i/i — r) + 7/iz] 



djdz 
Y^djdz 



(182) 



<c(e)/i 2 - £ - 55 / 2 0VT^(y-^)- 



Using inequality (h — r yz)h — ih = (n — i)h — jzh > h(l — jz) > h(l — 7) we obtain 
that 

|/J"| < Chh- b&/2 (j) VT ^{y -v) [ z 2 [ d 1 

Jo Jo 



E * 

{i-.T+h 5 /±<ih<T-h} 



(183) <Ch 1 - 5S/2 (f) VT=¥ (y-v). 

Now from (jEg) , (fTTKl) . (fTTHjl . (fT871 . (JTH2J) and (gEZD we obtain that 

(184) \D»(p®H-p ®' h H) (r, T, v,y)\< CW<\>^- T (y - «) 

for some positive 7. The last summand in the right hand side of (117411 admits the same 
estimate (I184JI because of the smoothing properties of the operation <g>' h . Hence, 



(185) 



\D v v {p - Pa)(t, T, v,y)\ < Ch^^(y - v) 
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Making a chaige of variables v = z + v' in ()170ft we get that the integral w.r.t. v is 
equal to 

cT 



dr / p(u, r,z,z + v')xu(r, v)D y v (p - p A ) (r, T, z + v', y)dv' 



(186) D v z 



Taking into account ()185|1 and making an integration by parts in (|186|) we obtain that 
(I186D does not exceed 



(187) 



ch? 



dr 



From (ITTOl) and (IT87T) we obtain that 



cw 

( 

\JT - u 



|/3|<^ 1+7 </yr(2/-^)- 

The estimate for Ii may be proved absolutely analogously to the estimate for I 3 . Thus, 
we proved that 

(188) hp <g> T x \p ® T x \p - p A }] = 0(/i 1+ Vvt(2/ - *)) 

The estimate 

h 1/2 p ®F x \p- p A ] = 0(/i 1+ Vvt(2/ " «)) 

may be proved by using the same decomposition of p — Pa- This completes the proof 
of Theorem 1. 
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